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A comprehensive computer code (FULLNOZ) has been developed to perform detailed 
calculations of rocket nozzle flows downstream of the sonic line. The code uses 
the streamtube method to integrate the hyperbolic governing equations of steady 
supersonic flow. The proyram represents a sipnificant advance in nozzle flow 
predictions through the coupling of gas/particle nonequilibrium effects, non- 
equlitbrium chemistry, turbulent boundary layers (to determine the effects of 
wall heat transfer and shear stress) and turbulent mixing across streamtubes. 
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This report describes the analytical and numerical techniques employed by the 
code, presents results of a sample calculation for the Minuteman Stage 2 nozzle 


and gives complete instructions on the preparation of input data and a iull 
FORTRAN listing. 
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SUMMARY 


This report describes the analytical and numerical techniqu..: utilized in 
the development of a fully-coupled rocket nozzle flow computer code (FULL- 
NOZ). The code uses the streamtube method to integrate the governing equa- 
tions of steady supersonic flow. The elliptic Navier-Stoxes equations are 
reduced to hyperbolic form by assuming diffusional effects along streamlines 
are small compired to those across streamlines, Finite difference techniques 
are then used to solve the hyperbolic equations along and perpendicular to 
streamlines. 
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FULLNOZ represents a signif.cant advance in nozzle flow calculations by 
coupling the effects of nonequilibrium chemistry, gas/particle thermal and 
dynamic nonequilibrium, turbulent mixing across streamtubes and turbulent 
boundary layers. Turbulent mixing is treated via a phenomenological eddy 
viscosity model, while the turbulent boundary layer analysis utilizes the experi- 
mental data of Keener and Hopkins (which relates the compressible skin fric- 
tion coefficient to measured velocityAemperature profiles in flows with favor- 
able pressure gradients), the Van Driest transformations and the momentum 
integral equation. The operation of FULLNOZ requires the specification of 
initial gas and particle properties just downstream of the sonic line, the nozzle 
wall contour and temperature, and a chemical reaction mechanism and rate 
coefficients. The marching scheme proceeds downstream computing flow 
properties and composition along surfaces orthogonal to a specified number of 
streamtubes, A mixed explicit/implicit differencing scheme is used to obtain 
the most favorable integration step size. 
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Sainple calculations are presented for the Minuteman Stage 2 nozzle. In 
addition, this report describes the preparation of input data and gives a full 
FORTRAN listing of the progra.n. .'41so0 included is an analysis of hete -oge- 
neous electron/ion recombination on particles, although it has not as yet been 
incorporated into the code. 
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PRE& as E 


This is a report on the work performed on Contract F44620-74-C -0006 
covering the period | September 1973 to 30 April 1975. It is basically a 
program user's manual for FULLNOZ, although some additional results ob- 
tained during this period, e.g. a preliminary study of how to incorporate 
heterogeneous electron-ion recombination into the code, have also been 
included. 


The authors would like to thank J.T. Kelly for his initial work in the 
development of FULLNOZ and Capt. L.R. Lawrence, AFOSR Program 
Manager, for his interest, encouragement and support during the development 
of FULLNOZ. 


This code is available for public use and may be obtained by forwarding 
a request and a tape to Lt. Robert Sperlein, DYSP, at the Air Force Rocket ' 
Propulsion Laboratory, Edwards Air Force Base, CA 93523, ‘ 


This scientific report has been reviewe. and is approved. 
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LIST OF SYMBOLS 


streamtube area; also pre-exponentiai term in rate coefficient 
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3 equations 
q 4i, jj enthalpy-temperature polynomial coefficients, see Eq. (20) 
B activation energy 
Ch particle drag coefficient : 
§ Cp normalized drag coefficient, see Eq. (34) : 
Copy _ free molecular drag coefficient 
| Cp, high Reynolds number drag coefficient, see Eq. (35) 
q Cj mass fraction of ith species 
P Cr skin friction coefficient, see Eqs. (54) and (59) 
’ Cp gas specific heat i 
Cz particle specific heat ; 
3 D diffusion flux, defined by Eq. (23) 
j D diffusion coefficient 
F; defined as C;/M; 


ratio of actual particle drag coefficient to the drag coefficient for 
Stokes flow 
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Bp ratio of actual particle heat transfer coefficient to heat transfer 
coefficient for Stokes flow 
H gas stagnation enthalpy 
Hi2 shape factor, defined by Eq. (69) 
He hat of formation at 298°K 
4 h gas static enthalpy 
h! particle/gas heat transfer coefficient 
4 I total number of gas species 
§ J total number of particle groups 
; Kp equilibrium constant 
L & k thermal (eddy) conductivity 
he forward reaction rate coefficient 
Kp molecular thermal conductivity in particle/gas interaction terms 
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total number of species on left or right side of reaction 


specific heat polynomial coefficients 
enthalpy constant of integration, see Eq. (31) 
entropy constant of integration, see Eq. (32) 
Cc Dp 

turbulent Lewis number, Le = ; 
Mach number 
molecular weight of ith species 
streamtube mass flow, defined by Eq. (13) 
number of reactions 

h’ r 


kg 


Nusselt number, Nu = 


SaaS 2 id 


distance normal to streamline, also exponent in Eqs. (57) and (58) 


turbulent Prandtl numbe:, Pr = ae 

laminar Prandtl] number 

static pressure 

heat flux, defined by Eq. (22) 

wall heat transfer, defined by Eq. (62) 

universa] gas constant, also value of nozzle radius at each x 
Reynolds number based on streamtube width 

particle Reynolds number 

Reynolds number based on distance along nozzle wall 
Reynolds number based on boundary layer momentum thickness 
radial distance from axis; also recovery factor, Eq. (49) 
particle radius 

distance along streamline 

entropy 

Stanton number, defined by Eq. (62) 

static temperature 


recovery temperature 


velocity along streamline 

velocity in boundary layer 

friction velocity, defined by Eq. (63) 
particle velocity normal to gas streamline 
production rate of ith species 

axial distance from nozzle starting line 
distance from nozzle wall 


molar concentration 


Greek Letters 


numerical stability coefficient for marching scheme, see Eq. (29) 
ratio of specific heats 

Gibbs free energy 

change in % across streamtube 

boundary layer thickness 

displacement thickness 

finite difference mesh spacing in n direction 

finite difference mesh spacing in s direction 


angle between streamline and plume axis; also boundary layer 
momentum thickness 


respectively; final, initial and intermediate streamline radii of 
curvature 


T/1000, used in Eq. (30) 

kinematic viscosity 

stoichiometric coefficient on left hand side of reaction 
stoichiometric coefficient on right hand side of reaction 
eddy viscosity 

molecular viscosity in particle/gas interaction terms 
gas density 

particle cloud density 

density of liquid or solid particle 


shear stress, defined by Eq. (21) and Fa. (60) 
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3 > summation | 
2 Subscripts 
aw adiabatic wall 
c compressible F 
e wall streamtube . 
7 FM free molecular 
g gas 
4 i ith species; also, incompressible 
j pavticle group identification index : 
e k streamtube index 
4 y orthogonal surface index; also refers to species on left or right ; 
hand side of reaction in Eq. (6) 
; m mth reaction 
, an differentiation in direction norma! to streamlines 
° stagnation value 
p particle 
8 differentiation in streamline direction 
t totai 
E w wall : 
4 
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I, INTRODUCTION 


This report describes a new fully-coupled rocket nozzle code (FULLNOZ) 
which treats, simultaneously, gas/particle nonequilibrium, nonequilibrium 
chemistry, diffusion across streamlines, and turbulent boundary layers within 
axisymmetric and two-dimensional nozzles, Because of its fully-coupled capa- 
bility FULLNOZ can be applied to a wider range of problems than such existing 
rocket nozzle codes,t as TDK! (which does not account for gas/particle non- 
equilibrium effects), the nozzle portion of the CONTAM® cede, the Original con- 
stant y gas/particle nonequilibrium code of Nickerson and Kliegel® and the 
recently-developed gas/particle and chemical nonequilibrium code developed at 
Lockheed/Huntsville.* The emphasis in the present code has been on an accu- 
rate calculation of nozzle exit plane gas and particle properties (particularly 
major and minor neutral and charged species concentrations), rather than on 
the determination of specific impulse, although FULLNOZ is well equipped to 
calculate Isp. 


t e.g. none of these codes calculate turbulent boundary layers. 


1. 'ICRPG Two-Dimensional Kinetic (TDK) Nozzle Analysis Computer Pro- 
gram,'' Dynamic Science Corp., December 1973 (revised version). 


2. Hoffman, R.J., English, W.D., Oeding, R.G., and Webber, W.T., 
"Plume Contamination Effects Prediction: The CONIAM Computer Pro- 
gram," Final Report, Air Force Rocket Propulsion Laboratory, AFRPL- 
TR-71-109, December 1971. 


3. Nickerson, G.R. and Kliegel, J.R., “Axisymmetric Two-Phase Perfect 
Gas Performance Program," TRW Systems Report No. 02874-6006-R000, 
Vols. land U, April 1967. 


4. Penny, M.M. and Smith, S.D., "Supersonic Gas-Particle Flows, Includ- 
ing Reacting Chemistry, '' JANNAF 8th Piume Technology Meeting, 
Colorado Springs, CO, July 1974. 


i 
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FULLNOZ is based on the MULTITUR3E code developed by Boynton, ? 

which incorporates the streamtube method! (described by Boynton and Thomson‘) 
to integrate the hyperbolic governing equaticns of steady supersonic flow. The 

4 major routines incorporated into FULLNOZ which are not in MULTITUBE in- 
clude (1) particle/gas nonequilibrium, (2) :onequilibrium chemistry and 
(3) turbulent wall boundary layers. Briefly, in the streamtube method the 
elliptic Navier-Stokes equations} are reduced to hyperbolic form by assuming 
that diffusional effects along streamlines are small compared to diffusion 
across streamlines. This assumption is very good for rocket nozzle (and plume) 
flows and enables one to solve an initial value problem (where a marching pro- 
cedure can be used) rather than the more difficult boundary value problem. The 
gas flow equations, in finite-difference form, are solved along and perpendicu- 
lar to streamlines while a full continuum particle cloud system of equations is 
incorporated for the condensed phase. 


sks 


The advantages of using the streamtube method over the method of charac- 
teristics in calculating rocket nozzle (and plume) flows are: 


e Species diffusion, shear, and heat transfer normal to streamlines 
are easily included. 


Chemical reactions or internal relaxations are easily incorporated 
since the calculation follows streamlines. 


* Bounding surfaces and gradients normal to streamlines are treated 
without difficulty; and 


¢ Awide variety of boundary conditions including mass transfer, shear 
and heat transfer can readily be incorporated. 


t In contrast to the abovementioned codes,! “4 which all use the method of 
characteristics. 


, Tn using the Na‘ ier-Stokes equations as a base the technique can readily 
. be extended, vith the proper boundary conditions, to low density nozzle 
: flows. 


5. Boynton, F.P., "The MULTITUBE Supersonic Flow Computer Code," 
General Dynamics/Convair GDC-DBB 67-003, February 1967. 


6. Boynton, F.P. and Thomson, A., "Numerical Computation of Steady, 
Supersonic, Two-Dimensional Gas Flow in Natural Coordinates," J. 
Computational Phys. 3, 379-398 (1969). 
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Figure 1. FULLNOZ schematic. 


Operation of the code is achieved by specifying (1) initial gas and particle 
properties in the supersonic region just downstream of the nozzle throat, 
(2) the nozzle wall contour, (3) a chemical reaction mechanism and rate coeffi- 
cients, (4) physical properties of the particles, and (5) the nozzle wall tempera- 
ture (for the boundary layer calculations), With the above input data the march- 
ing scheme steps from one orthogonal surface to the next (Fig. 1) computing 
gas and particle properties within a specified number of streamtubes. The 
(essentially) continuous particle size distribution is represented by up toa 
maximum of 8 discreet particle sizes. All particles cross the gas stream- 
lines, but the lighter particles follow the gas streamlines more closely than the 
heavy particles. As depicted in Fig. 1, limiting particle streamlines are com- 
puted for each particle size. At each orthogonal surface (i.e. each integration 


It would be useful to incorporate an "initializing" scheme which utilizes 
combustion chamber properties as initial conditions, computes flow 
properties through the transonic region and estabtishes an initial super- 
sonic data line. Such a scheme, developed by Nickerson and Kliegel,? 
is contained in the codes of Refs. 1-3, and with additional programming 
could be inserted into FULLNOZ. 


1] 


E step) the code calculates wall shear stress and heat transfer, boundary layer 

¥ displacement thickness and velocity and temperature profiles. The shear stress 

3 and heat transfer are coupled to the main nozzle flow via their effect on the 

a wall streamtube properties.! The turbulent boundary layer analysis utilizes 

E the Van Driest’ transformations and the experimental data of Keener and 
Hopkins, *® which relates the compressible skin friction coefficient to neasured 
velocity/temperature profiles in flows with favorable pressure gradients (see 

g Section I1.C.4). The boundary layer momentum thickness is computed via the 

E momentum integral equation, using the wall streamtube pressure, velocity, etc., 

3 as the boundary layer ‘edge! conditions. The displacement thickness is then 

determined from the velocity profile and momentum thickness. 


The code will not handle shocks that might originate from the nozzle wall 
in turning the flow. The flow in the region of the shock will be treated as a 
strong compression wave. i 


t This implicitly assumes that the boundary layer momentum and energy 
thicknesses are smaller than the wall streamtube thickness. This will 
generally be true for high Reynolds number nozzle flows. Provision is 
in the program, however, to transfer momentum, heat and mass across 
streamtubes, so that if an appropriate transfer coefficient can be defined 
the boundary layer effects can be "felt" throughout the flow. 


Nozzle shocks could be incorporated into FULLNOZ in a manner similar 
| to that employed in the AIPP code?! to detect internal plume shocks, but 
this would require additional programming. 


«J 
e 


Van Driest, E.R., "Turbulent Boundary Layer in Compressible Fluids," 
J. Aeron. Sci. 18, 145-160 (1952). 


oo 


Keener, E.R. and Hopkins, E.J., "Van Driest Generalization Applied to 
Turbulent Skin Friction and Velocity Profiles Measured on the Wall of a 
Mach 7.4 Wind Tunnel," AIAA J. 11, 1784-1785 (1973). 


9. Pergament, H.S. and Kelly, J.T., "A Fully-Coupled Underexpanded 
Rocket Plume Program. (The AIPP Code). Part I. Analytical and 
Numerical Techniques, '' AeroChem TP-302a, AFPRL-TR-74-59, 
November 1974, 

10, Pergament, H.S. and Kelly, J.T., "A Fully-Coupled Underexpanded 
Rocket Plume Program (The AIPP Code). Part U. Program User's 
Manual, "AeroChemTP- (in preparation). 
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The code was written for operation on a CDC 6400 computer and requires 
approximately 122 K (octal) of core storage. For operation on machines with 
larger core storage lhe user may want to increase the dimension of some param- 
eters, such as total number of species (25 maximum) and total number of 
reactions (40 maximum). 


This report serves primarily as a program user's aes but also con- 
tains all the governing equations (Section I) and the results of some preliminary 
parametric calculations for the Minuteman, second stage nozzle (Section IV. A), 
The user should be cautioned that, although the code has been formally debugged 
for several test cases, extensive calculations have not (as of April 1975) been 
made; therefore some operational problems may be experienced. If these do 
occur, please contact the authors. 


ll. GOVERNING EQUATIONS 


This section gives the governing differential ard finite-difference equa- 
tions used in the code. Also included are the auxiliary equations for calculat- 
ing thermodynamic and chemical kinetic properties of the system, zas/ 
particle drag and heat transfer coefficients and turbulent boundary layer 
properties. 


A. Differential Equations 


1, Gas Phase 


For most high Reynolds number nozzle flows of interest the turbulent 
transport of mass, momentum and energy throughout the main nozzle flow 
will be negligible. Thus in the equations that follow all transport terms will 
be identically zero. In practice this is achieved by setting ITURB (Card ¢, 
Cols. 56-60) equal to zero. If turbulent diffusion across streamtubes is to be 
included due either to initial non-uniformities or to the propagation of boundary 
layer effects across the flow, appropriate (constant) values of the eddy trans - 
port terms, up, Pr and Le will have to be input on Card 10. Ifa suitable expres- 
sion for edly viscosity in terms of local properties can be developed for trans - 
port across streamtubes within nozzles, this expression can easily be added to 
the program in subroutine TRANSP. 


t 


The same nurnerical integration techniques and similar program logic is 
; :0 ‘ ; . 
incorporated into the AIPP code," to which the reader is referred for 
additional information on program subroutines, 
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Global Continuity 


(pu), + ou HPP 4 pus, = 0 (1) 
s-Momentum 
J p f 
1 7 9 P; P. g ( 
puug + Py 5 Lf wunh-) 3 Por, * «mp ) (2) 
jst Pi 
n-Momentum 
J Pe f, 
9 P 
pve,+e,+) 3 —i-i v, (3) 
jet : of ; 
Species Continuity? 
BO vt aaa 1 Le 
palCi)s = pW: + a [= + 5, (Ci)n Ie (4) 


1 
J Ke . 
Bj j 2 
+) 7 oy M+ (uu ) 
2 Psp wee P 
j=: j 
Cc 
a. : (1 )! (5) 
= 7 . 
3 Pry P; { 


t 


In this formulation, Wj is expressed in units of sec™!, 


The drag (fp;) and heat transfer (gp;) factors are defined as, fp; = CpRe,/ 24 
and Bp; = (Nu/2) (Tp 96 T,W/(Tp - T), where T,, is the recovery temperature 


based on the relative velocity between the gas and particle. 


The following auxiliary expressions are required: 


a 


Species Production 


he sa ~ = i ts 
Ww; * ry ») (vr - vw) [Kien Pi Yy 
m1 is 
k L Y 


Equation of State 


(7) 
Cc; 
RT Poet Be 
M; 
eye 
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Stagnation Enthalpy 


H = = + h(T) (8) 


2. Particles 


For the condensed phases present within the flow a continuum parti- 
cle cloud assumption is made and therefore field conservation equations for 
continuity, momentum and energy can be written for the particles. The 
(essentially) continuous distribution of particle sizes at each point in the flow 
is modeled by several groups of constant size particles representative of the 
distribution. Fora given group, j, the conservation equations, written ina 
streamline oriented coordinate system are: 


Continuity 


Re me 


4 (x ), + (re, Vp Jy t+ FP 8 -rp VV 8.29 (9) 
‘* ba s 

J "2°? ae Poms BT Ye, Py 

3 s-Momentum 

P;"p, os + Vo _(u ) -V. (Vv. 6. -u @,) = 

4 pe ey 

on 

x 9 ( 

: ; ; (10) 

2 es rT, * "P; ) 
J 
n-Momentum 
=) = 
Wp (Vv, )s + Vp (Vp ), FY. (Ve n*%, 95) 
jj jj y 4 J 

q Bt 

‘ see ey (11) 
‘ 2 P.r. 2 

s “p, j 
J 
4 Energy! 
, “Ka Ep, fi, 
up (Cc, Ty ). + Vp. (Cs T) qe = ar? Sie ce - T) 
j j j J Pj s j 
t 


The effects of chemical reactions on the surface of particles are not 
included in this analysis. 


When the particle undergoes a phase change (liquid to solid) it is kept 
at the solidification temperature until the total heat of solidification is released 
(via radiative anc convective heat transfer) to the gas. 


B. Finite-Difference Equations 


A finite difference formulation of the gas phase and particle cloud govern- | 
ing equations is utilized on a grid which lies along and perpendicular to the | 
streamlines. 

| 
i 


1. Gas-Phase 


The momentum and energy equations (Eqs. (2), (3) and (5))are solved 
via an explicit finite-difference marching technique, whereas the species 
continuity equation (Eq. (4)), utilizes an implicit finite difference formulation 
developed in an earlier study at AeroChem."! This mixed form of the differ- 
ence equations is necessary for the economical operation of the present code 
since, for near-equilibrium chemistry, the explicit finite-difference form of 
Eq. (4) leads to stability-limited (impractically small) integration step sizes. 
Equations (2), (3) and (5) are left in explicit form since the required integra- 
tion step sizes for stability are reasonable; an implicit formulation of these 
equations would unnecessarily complicate the calculations. 


To minimize the effects of large tube-to-tube property variations on 
the calculation of streamline curvature, a Von Mises -type transformation is 
employed, i.e. 


in = an J pur dn (13) 


Using Eq. (13) the finite difference forms of the differential equations become: 


Global Continuity 


is a oe Eee ee 41 (14) 


ll. Mikatarian, R.R., Kau, C.J., and Pergament, H.S., 'A Fast Computer 
Program for Nonequilibrium Rocket Plume Predictions," Final Report, 
AeroChem TP-282, AFRPL-TR-72-94, August 1972. 
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The species mass fractions at station k,f i 1 (C; k, £+1) are determined 


by linearizinyg the chemistry terms ( (Wik) and inverting the resulting 
matrix (see Ref. 11). 


Energy 
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a. Particles 
The finite difference forms of Eqs. (9-12) are: 
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3. Integration Step Size 


The integration step size must be limited in order to perform a 
stable calculation. The stability of the explicit finite difference scheme for 
solving the gas dynamic equations is discussed by Boynton and Thomson,°® 
who show that the stable step size for laminar flow is determined from, 


1 1 (28) 
ti ? [ de i 2 i/ | 
(M? - 1)°? 


For turbulent flow the eddy viscosity replaces the laminar viscosity in the 
expression for Re.! The step size determined from Eq. (28) may not be 
sufficiently small if the chemistry is very "fast". Thus, an additional 4 
"chemistry" control has been incorporated into the code: at each integration ‘ 
step (orthogonal surface) the species mass fractions are tested for sign. If { 
any Cj goes negative the integration step size is halved and the calculations ¥ 
are repeated until the Cj in question becomes positive or the step size be- 

comes less than the minimum aliowable step size. In the latter case the j 
program terminates. { 


The particle cloud system of equations (Eqs. (9-12)) contains deriva- 
tives related only to the convection of mass, momentum and energy. In the ! 
momentum and energy equations, the convection terms are equal to the 
particle/gas interaction terms, Unlike the gas flow equations the particle 
equations have no wave or diffusive nature. Consequently when these equa- 
tions are solved explicitly there is no stability limitation on the integration 
step size and the step size can be determined from Eq. (28). 3 


4. Solution Procedure 


The gas phase and particle cloud finite difference equations are 
solved on a grid consisting of the streamtubes and the surfaces orthogonal to 
them, as illustrated in Fig. 1. All gas and particle flow properties, stream- 
line positions and angles must be known along the initial orthogonal surface 
downstream of the throat (M >1). Starting with the first streamtube (adjacent 
to the axis) the streamtubes are extended a distance 6s to a downstream sur- 
face, utilizing the radius of curvature, (Qs)! obtained from the normal momen- 
tum equation (Eq. (16)) and the known initial pressure distribution along the 


t 


Eq. (28) is not, of course, valid in the limit, Re 0. In that case (i.e. 
2 


rs) i 
the situation for most practical nozzle flows), 6s < > fM?-1} ‘ 


surface. The resulting streamtube areas are then used (o determine all the 
necessary gas propertics at the downstream surface from Eqs. (14), (15), 
and (17)-(23). The transfer of mass, momentum and energy into or out of 
each streamtube is arranged such that what is lost froma given streamtube 
is gained by the adjacent tube. Thereby mass, momentum and energy are 
automatically conserved. 


In order to render this marching scheme conditionally stable a single 
iteration on the radius of curvature calculation is required,° with the new value, 
(8,)", determined using the downstream surface properties. (9©,)” is then 
combined with (9,)’, in the following expression, 


(8,) = (1 - @)(@s)’ + a(85)” (29) 


This new value for the radius of curvature is then used in the calculation. A 
value for a of 0.55 is used in the present code. (The scheme is conditionally 
stable for a 2 1/2). 


After the gas phase equations have been solved for the first stream- 
tube the particle properties are determined by sequentially applying the parti- 
cle momentum, continuity and energy equations. The calculation of gas and 
particle properties then moves outward to the next streamtube and the proce- 
dure is repeated up to the last streamtube! where boundary conditions must 
be aprlied. 


C. Auxiliary Calculations 


1. Thermodynamic Input Data 


The thermodynamic data are input via curve fits of specific heats!® 
of individual species in the JANNAF tables.!? These curve fits have the form, 


Cy = Ly t+ Le ® + Ls. 07 414 G2 45, 07? cal/mole-°K (30) 
1 


t In general the limiting particle streamJines for each particle size will 
not extend to the last streamtube,. 


12. Cruise, D.R., "Information Manual for the Theoretical Propellant 
Evaluation Program," Naval Weapons Center PEP NOTE TN-U-1 (plus 
additions), December 1964. 


13. JANNAF Thermochemical Tables (Dow Chemical Company, Midland, 
Mich.), continuously updated. 
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T 
h = { ee dT + Ly, kcaY mole (31) 


"7298 
where L,. = H? - C, dT; H® being the heat of formation at 298°K. 
i f Pp f 8 
0 i 


The entropy is expressed as, 


: dT fe) 
S, = { Cp rT + Ly, cal/mole--K (32) 


where L;; is the entropy integration constant. The coefficients L, - L; are 
input on Card group 8. 


2. Chemical Kinetic Input Data 


Ten possible reaction types are included in the program: 


a where ® = (°K)/1000. The enthalpy is then expressed as, 
Reaction Type 


(1) A +B z C+D 
(2) A+B+MzC+#t+M 
(3) A+B z= C+tD#teE 
(4) A+B = uC 
(5) A+M z= C+tDtM 
(6) A +bB —- ¢c +obD 
4 (7) A+t+Bt+M-mrC+#+M 
4 (8) A+B -~C+DHE 
(9) A +B = ¢ 
A+M -C+DtM 


| (10) 


Reaction types (6)-(10) correspond to reaction types (1)-(5), but proceed in 
the forward direction only, In Reactions (2), (5), (7) and (10), Mis an 
arbitrary third body. In this program, all species are assumed to have equal 
third body efficiencies, 


2+ 


SIFY SET LE GT EE TOT Te PITY STE ewe FPR aR rave ye ee ee ra ee i teiccenae Mi Beat otal ai tla 3 


ee sceneries a te ee eee 


The forward rate coefficient, ks, is input to the code as one of the 


following 8 types 


Rate Coefficient Type 


(1) 
(2) 
(3) 
(4) 
(5) 
(6) 
(7) 


(8) 


cor Sere eo 
" 


The equilibrium constant, Kp) 


In K, 2 


A 
ar” 
AT” 
1 
AT-2 
A exp(B/RT) 
AT exp(B/RT) 
AT 2 


atN exp(B/RT) 


is determined from 


- AG/RT (33) 


where the Gibbs free energy, AG, for individual reactions is computed from 


the input thermodynamic data. 


T Rate coefficient data for typical rocket nozzle and plume reactions may 
be found, e.g., in Ref. 14. 


14. Jensen, D.E. and Jones, G.A., "Gas-Phase Reaction Rate Coefficients 
for Rocketry Applications,'' Rocket Propulsion Eotiblishment Technical 
Report No. 71/9, October 1971. 
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3. Particle/Gas Drag and Heat Transfer Coefficients 


The momentum and energy exchange (via convective heat transfer) 
between the small diame‘er particles and the combv stion products 1n the 
nozzle cannot adequately be described by simple theoretical expressions (e.g. 
Stokes law).!5 Empirical correlations of drag and heat transfer coefficients 
developed by Crowe!® have therefore been incorporated in the particle/gas 
interaction terms (see, e.g. Eqs. (2), (3), and (5)). 


a. Drag Coefficient - The drag coefficient has been correlated by Crowe’® 
in terms of a normalized value, 


Cp = (Cp - Cp V/(Cdpyy - Spy) (34) 


where Cp, is the drag coefficient at very large Reynolds number and Corm is 


the free molecular drag coefficient. Note that for Rep << il, Cp ~ 1, while 
for Rep >> 1, Cy 0. 


The expressions needed to evaluate Cp from Eq. (32) are: 


Cp, = 0.66 + 0.26 [exp(4 InMp)-1] + 0.17 exp[ -2.5(InMp/1.4)*] (35) 


_g2 4 2) 
| Cy . exp(- 5,°/2) (14 25,2) + HS-+8 pet erf(S,) + < fi (36) 
FM Nr S) 2s,* N T/T, 
. where, S; = Vy/2 Mp (37) 
Cp = G(Kn) D(Kn, Rep) (38) 


0-4 0.5 
G(Kn) = Kn?" exp(1.2 Kn’ "”) 


(39) 
1 + Kn°* exp(1.2 Kn") 


15. Soo, S.L., Fluid Dynamics of Multiphase Systems (Blaisdell Publ. Co, 
Waltham, Mass., 1967). 


16. Crowe, C.T., "On the Momentum and Heat Transfer Equations for Two- 
Phase Plumes,"' Washington State Univ,, March 1971. 
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24/Re,, (42) 


Qa 
0 
° 
uw 


The expression used in the particle/gas interaction terms is then, 


Cc 
D 
f. = FFF 
P Cy, 


(43) 


where FFF is a factor (input on Card 15, Cols. 1-10) used to arbitrarily vary 
Gp to account for uncertainties in the above analysis. 


b. Heat Transfer Coefficient - Heat transfer from the particle to the 
gas is expressed in terms of a Nusselt number as, 


q = awry Nu k(T, - T,) (44) 
where T,, is the recovery temperature. From Crowe!® we get the following 
expression, 

2 Be -Re 
Nu = Nuyp + - Re,Pr, exp( Re,/2Mp) (45) 
sy 
Nugp = Nuo/( 1 + SPT (Mp/RepPr,) Nug ) (46) 


Nu, is the Nusselt number in incompressible flow, expressed as 


Nu = 2.0 + 0.459 Rep’*? Pr*3 (47) 
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The recovery temperature is defined as, 


Tr = T + r( (u-up? + V2 )/2Cp (48) | 


where the recovery factor r is!®; | 


r= 0.9 + (rey, - 0.9) exp(-Re,/2M,) (49) 
and 
tem? a (2 +0. 67 exp(-M‘4) ) (50) 


The expression used in the particle/gas interaction term is, 


Nu s me 
g, = FFG ToT (51) 


where FFG is a factor (input on Card 15, Cols. 11-20) used to arbitrarily 
vary 8p to account for uncertainties in the above analysis. 


D. Turbulent Boundary Layer Equations! 


The turbulent compressible boundary layer analysis is initiated by 
calculating the corresponding adiabatic flat plate, zero pressure gradient, 
incompressible boundary layer properties, based ona 1/7 power law!” 


1 


The "free stream" properties (subscript "e'') for the boundary layer 
analysis are taken to be the properties in the last (wall) streamtu e. It 
is implicitly assumed that the boundary layer displacement thickness is 
smaller than the width of the wall streamtubc. 


17. Schlicting, H., Boundary Layer Theory, 6th Ed. (McGraw-Hill, New 
York, 1968), p. 599. 


Were ses ener ee 


sg ee Te een 


velocity profile.! Expressions for the incompressible momentum and bound- 
ary layer thickness are, 


8 s(x) 0.036/Re,°* (52) 


6:(x) = 10.286 @ ;(x) (53) 


The incompressible skin friction coefficient is evaluated from the Karman- 
Schoenherr relation. *° 


Cy. = [17.08 (logio Reg. + 25.11 logio Reg, + 6.012]~ (54) 


Transformation from the compressible to the incompressible boundary 
layer is accomplished via the Van Driest equations,’ which relate the velocity 
and temperature profiles as follows, : 

T -T 
= e w U 
oe ee ‘( a (55) 
T u.) 


T it was originally anticipated that the complete boundary layer equations 
including nonequilibrium chemistry would be solved by finite differences 
and coupled to the nozzle flow solution. A number of boundary layer 
codes are available, including tnose of Herring and Mellor!® and the 
Aerotherm BLIMP code,!? which would be very useful for this purpose, 
However, initial attempts to incorporate the BLIMP code into FULLNOZ 
showed that it would take more effort than was warranted at this time. 
Consequently, a more simplified analysis was incorporated into the 
present code. 


18, Herring, H.J. and Mellor, G.L., "A Method of Calculating Compressible 
Turbulent Boundary Layers,'' NASA CR-1144, September 1968. 


19. Tong, H., Buckingham, A.C., and Morse, H.L., "Nonequilibrium 
Chemistry Boundary Layer Integral Matrix Procedure," Aerotherm 
Final Report No. 73-67, July 1973. 


20. Hopkins, E.J., Keener, E.R., and Louie, P.T., "Direct Measurements 
of Turbulent Skin Friction on a Nonadiabatic Flat Plate at Mach Number 
6.5 and Comparisons with Eight Theories, '' NASA TN D-5675, February 
1970. 
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RTE NT I EEL TPT a i ne ee eee = wee: 


Fe 


| T, = (= ) (56) 


Equation (56) becomes the Crocco relation for n = 1 and the quadratic for 
ns 2. 


A general expression between T/ T, and u/ U, can then be written 


U 
: Rea i BCG ys (57) 
where 
“se *, a T, 
As(=2-iibe (et -P)ices 


From Back and Cuffel,”! wall cooling (Ty < Tay) tends to make the 
exponent n closer to 1.0, while wall heating causes n to approach 2 or more. 
The value ofn selected here was 1.2. It was determined by fitting recent 
data of Keener and Hopkins® for which, Te) Taw = 0.32. (This corresponds 
to a wall temperature of 1090-1300 K.) Equation (57) was used to generate a 
table of T/ Te Vs. U/Ue, which is then used to determine the compressible i 
skin friction coefficient via the Van Driest transformation 


Cc 
f. 2 
e, : ‘ e :) a( U, } (58) 
The Stanton number was determined fror: the relation 
S; = 0.35 te (59) 


which is supported by the experimental data of Back and Cuffel* for accelera- 
ting flows. 


21. Back, L.H. and Cuffel, R.F., "Relationship Betwcen Temperature and 
Velocity Profiles in a Turbulent Boundary Layer alors a Supersonic 
Nozzle with Heat Transfer,'' AIAA J. 8, 2066-2069 (1970), 


; 


Wall shear stress and heat transfer are then calculated from 


1 2 
Tw = 7 Pe Ue ee (60) 
Gy = ~St Pe Ue Cy (Ty, - Ty) (61) 


In order to obtain the velocity profiles the friction velocity profile, 
(U/U +); = {(U/U,) is first determined from another Van Driest transformation,’ 


/, u/U, Ts 
(Exe (ey L" (BYE) ow 


where the friction velocity is defined as, 
3 Y, 
Ur = (be Ue Cy/2py, ) (63) 


(yf); is then obtained from the standard incompressible boundary layer 
profiles 


T U U 
ge Ur, ; wee 5 (64) 
yU U/U_ + 3.05, 
— = exp gael Iie ; 5 8 ~ = 13.96 (65) 
( : 
yU (u/u,- 5.05) 
artis exp as H Ue 13.96 (66) 


' The values of yU+/v are normalized by exp{[U/Ur7 max 75-05]/2.5} to 
avoid the calculation of the v profile in the boundary layer. 


22. Kays, W.M., Convective Heat and Mass Transfer (McGraw-Hill, New 
York, 1966). 
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Another Van Driest transformation yields (y/5). = (y/6);. Once the com- 
pressible profiles are known, @/6, 5*/5, and the shape factor are determined 


from 


e/6 = { (a )0- a) ah) (67) 
w= S, (a) a4) s) 


6*/6 
Hn = 675 (69) 


The variation of momentum thickness along the nozzle is determined via 
: integration of the momentum integral equation, 


Cc 
f du ap 
| ¢ , 4 1 Me, 9 Me, o aR 
5 = 3, + [Het2je U; ax ars a oe ae (70) 


3 by the use of backward differences.’ 6 and 5* are then evaluated from Eqs. 
(67) and (68). 


| Il, PREPARATION OF INPUT DATA 


All necessary information for preparing input data is given below; Fig. 2 
defines some of the input for a sample case. Many of the input parameters 
that are left blank are used in the companion rocket plume code!® (AIPP), but 
not in FULLNOZ. 


t This requires that a throat value of O be assumed. 


ne emelebstnt hea, Ste atk - 


| 
i 
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Z Ly x8(3) 
te yyyyyra, x18 x8) x82) RB(3) 
LZ Z-W8\+RB() _RBA2) 
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KMAX®8 
NBL(J)*7 ONOSP (u) 
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CONTAINING PARTICLES 
OF SIZE GROUP J 
LIMIT OF 
PARTICLE CLOUD 
OF SIZE GROUP J 
TREDY (u) 


Figure 2. Input data definitions. 


Card Fortran 
No. Cols. Name 
1 1-10 NDATA 
2 1-80 ID 
3 1-10 NSEC 
11-22 DXLSS 


Description 


Number of data sets; in present ver- 
sion of code data sets cannot be 
stacked; therefore NDATA = 1 


Run identification 


Maximum run time (sec); when run 
time reaches NSEC nozzle properties 
at last orthogonal surface are punched 
and can be used to continue the calcu- 
lation, (Not operational in present 
version of code ) 


Approximate distance (cm) along axis 
between orthogonal surface print sta- 
tions. For control of print increment 
by number of integration steps (KP, 
Card 5, Cols. 1!-15) instead of axial 
distance, make DXLSS larver than 
XLMAX (below) 


Format 


I. 10 


20A4 
I. 10 


El2.5 


Fn nee Re 


Card 
No. Cols. 


q 3 23-34 


11-15 
16-20 
21-25 


26-30 


31-35 
36-40 
41-45 


46-50 


51-55 
56-60 


eS SLI ARON oo sen 


Fortran 
Name 


XLMAX 


ITYPE 


IKIND 


IBUGSH 


IPART 


ITURB 


KMAX 


NN 


Description 


Maximum distance (cm) along axis for 
which calculations will be made. 
XLMAX should be set somewhat larger 
than axial distance to nozzle exit plane 
in order to complete calculation at last 
wall point 


Inner boundary condition (nozzle axis) 
Set ITYPE = 1 


Outer boundary conditicn; set IKIND=1 
for wall boundary condition 


Not used; leave blank 
Not used; leave blank 
Not used; leave blank 


Debug printout index; 
0 - No debug printout 
1 - Extensive printout for debugging 


Not used; leave blank 
Not used; leave blank 
Not used; leave blank 


Particle indicator; 
0 - No particles 
1 - Particles in flow 


Not used; leave blank 


Turbulent flow indicator; 

0 - Inviscid 

3 - Constant values of turbulent mix- 
ing parameters, »p, Pr and Le are 
input on Card 10 


Initial number of streamtubes plus one 


(axis is counted as K=1; maximum is 
40) 


Number of terms for Cp polynomial 
curve fit plus one; for curve fits sup- 
plied with program,'? NN=6, If fewer 
coefficients are used for a given 
species NN is Icit at 6 and zeroes are 
input for the missing coefficients (see 
Card 8) 


34 


Pormat 


E12.5 


15 


15 


15 


I5 


te eA hunts one 


a te re Aa ED 


NR OE ET A Ee EE EE 


Bt Dagon na. catznston aN ES eB ote SL = PTT 
A Pe is wo Sal ae ss aes. cs 3 bos Bie Ss A 


eo Seen aoe S's ae oe cece enuten neat erent ened tne ft EAPO LIN ADAPT TAA 
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“i Card Fortran 

g _No. Cols. Name Description _ Format 

4 

e 5 11-15 KP Output control; number of orthogonal : 
3 surfaces between print stations, To 

3 control print via axial distance (DXLSS, 


Card 3, Cols. 11-22) instead of orthog- 
onal surfaces set KP larger than 
LPLANE (Card 5, Cols. 36-40) 


16-20 MMAX Number of points needed to describe 5 
shape of inner boundary; always set 
equal to 2 for axis boundary 


21-25 NMAX Number of points required to define 15 
nozzle wall contour (50 maximum) 


26-30 NDS Total number of gas species in flow 15 
(25 maximum) 


31-35 NITER Maximum number of iterations allowed 15 
for iterative solutions (e.g. calculation 
of streamtube propertics for variable 
Y) 
Recommended value: NITER = 50 


| 
| 
| 
| 


q 36-40 LPLANE Maximum number of integration steps [5 
3 for entire calculation 


41-45 IKINE Number of chemical reactions tS 
(40 maximum) 


NOTE: IF NO BOUNDARY LAYER CALCULATIONS ARE TO BE 


§ MADE CARD 6 MAY BE LEFT BLANK, BUT MUST BE 
INCLUDED : 
i 
6 1-5 TWALL Nozzle wall temperature (°K) £5.0 ; 
(assumed constant) : 
8 IBLFLG Boundary layer property printcut I} 
indicator; i 
0 - Print Re, 6,, Oj, Ce, Uq, 4, 8, 


Oo Cre: Ho, Aw ere 
1 - Print above plus velocity and tem- 
perature profiles 
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Card Fortran 
No. Cols. Name _ 
6 1l IBL 
indicator; 


AND NDS SPECIES 


8.1.1 1-13 A(L) Li 
14-26 A(2 L 


40-52 A(4 1 


obtained for flows without particles. 


r 36 


Boundary layer calculation 


0 - Do not calculate boundary layer 
1 - Calculate boundary layer 


7 1- 9 ALPHAH Factor that multiplies maximum stable 
step size (see Eq. (28)); Recommended 
value: ALPHAH = 0.8 


10-18 EPSLON Amount by which streamtube Mach num- F9.4 
bers must exceed one in order for the 
calculation to continue. 
value: EPSLON = 0.01 


Recommended 


19-27 TOL Convergence tolerance for iterative 
solutions; Recommended value: 
TOL = 1x10" 


28-36 DELTA Metric exponentt; 
0 - Two-dimensional flow 
1 - Axially symmetric flow 


37-45 ATOL Maximum allowable fractional change 
in streamtube area per step. Recom- 
mended value: 


) Specific heat polynomial con- 
‘ tants for first species, see 
27-39 A(3 Ls: 2 
( > Section I1.C.1 (cal/mole °K) 


t The particle conservation equations are written only for axially 
symunetric flow. Thus two-dimensional solutions can only be 


Format 


Il 


F9.4 


F9.4 


F9.4 


F9.4 


NOTE: CARD GROUP 8 DEFINES THE THERMODYNAMIC DATA 
FOR EACH SPECIES. THERE ARE 2 CARDS PER SPECIES 


E13.5 


Card Fortran 


No. Cols. Name Description Format 


8.1.2 1-13. A(5) Ls, Specific heat polynomial con- E13.5 
stant for first species 
14-26 A(6) Ly, Enthalpy constant of integration E13.5 


for first species 


° 298 
( AH, - ‘ Cy aT ) (kcal/mole) 


£6 eT OR OE Ee NED wth IEE esate RAT al tat 


oe ee eS 


27-39 = CS(1) L7, entropy constant of integration 
for firat species (cal/mole °K) 


8. NDS, 1 Repeat thermodynamic data for NDS 
8& NDS.2 species 


NOTE: CARD GROUP 9 IDENTIFIES EACH SPECIES AND 
SPECIFIES VARIOUS TRANSPORT PROPERTIES. 
MUST BE IN SAME ORDER AS CARD GROUP 8 


9.1 1- 4 IDENT(1) Species name; first species (the A4 
remaining data on Card 9.1 also apply 
to the first species) 
13-24 MUO(1) Viscosity at reference temperatuxe t E12.4 
(g/cm-sec); only used when flow 
contains particles, FOR IPART =0, 
MUO MUST = 0 FOR ALL SPECIES 
25-36 TO(1) Reference temperature for viscosity El2.4 
(°K) 
i 37-48 OMEGA(1) Exponent describing viscosity/tem- E12.4 
; perature relation; p « T 
$ 
49-60 PR(1L) Reciprocal of species Prandtl number E12.4 
61-72 SC(1) Reciprocal of species Schrnidt number E12.4 
73-80 MW(t) Species molecular weight E8.2 


t Values for common gases can be found in most physics and chemistry 
handbooks, 


Card Fortran 
No. Cols. Name Description Format 
5 
9.NDS Repeat for each species 


NOTE: CARD 10 IS NOT NEEDED IF ITURB = 0 (CARD 4, 
COLS. 56-60) 


10 1-10 TLE Turbulent Lewis number (constant) 
11-26 TPR Turbulent Prandtl number (constant) ; 
21-30 EDDYK Eddy viscosity, g/em-sec (constant) 


NOTE: CARD GROUPS 11 AND 12 DEFINE THE LOCATION OF 
THE INITIAL ORTHOGONAL SURFACE, AND THE FLOW 
ANGLE, PRESSURE, TEMPERATURE, VELOCITY AND 
SPECIES MASS FRACTIONS WITHIN EACH STREAMTUBE 


11 1-12 X(1) Initial axial position on inner boundary E12.4 
(cm) 
13-24 R(1) Initial radial position on inner bound- E12.4 
ary (cm); usually equal to 0 for axis 
boundary 
25-36 PHI(1) Initial flow angle! on inner boundary E12.4 
(radians); usually equal to 0 for axis 
boundary 
2m) I | 1-12 (2) Axial position at outer boundary of first E12.4 
| streamtube (cm) 
| 13-24 R(2) Radial position at outer boundary of El2.4 
first streamtube (cm) 
| 25-36 PHI(2) Flow angle at outer boundary of first El2.4 
streamtube (radians) 
37-48  P(2) Average pressure in first streamtube E12.4 
(atm) 
49-60 T(2) Average temperature in first stream- E12.4 
tube (°K) 
61-72 U(2) Average velocity in first streamtube Ei2.4 
(cm/sec) 
12.1.2 1-10 C(1,2) Average mass fraction of first species E10.3 


in first strearntube 


T Plow angle is detined as the angle between the flow velocity vector 
and the axis. 
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Card Fortran 
No. Cols. Name Description Format 


12.1.2 11-20 (2,2) Average mass fraction of second species E10,3 


in first streamtube ° 


71-80 C(8,2) Average mass fraction of eighth species E10.3 
in first streamtube 


CT ee SS Tae OT 


121.1828) 
NOTE: 12°) DENOTES THE NEXT INTEGER LARGER 
THAN (APS), REPEAT CARD 12.1.2 TO INCLUDE 
NDS SPECIES, 8 SPECIES/CARD | 
12.2.1 1-12 = X(3) Axial position at outer boundary of El2.4 
: second streamtube (cm) 
e e 
61-72 U(3) Average velocity in second streamtube E12.4 
(cm/sec) 
12.2.2 1-10 C(1,3) E10.3 
e e 
e e 
71-80 (8,3) : 
NDS ’ ‘ 7 
2.2. \(->- ) C(NDS,3) Average mass fraction of last species E10.3 
° in second streamtube 


e 
12. (KMAX- 1),1 1-12 X(KMAX) Axial position at outer boundary of last 612.4 
(KMAX-1) streamtube (cm) 


12. (KMAX- 1). oe ) C(NDS, Average mass fraction of last species E10.3 


: KMAX) in last streamtube 
NOTE: CARD 13 DEFINES THE POSITION OF THE INNER 
BOUNDARY 
13,1 1-12) XW(1) Initial axial position (cm); repeat of Bl2.4 
X(1) on Card 10 


13-24 RW(1) Initial radial position (cm); repeat of E12.4 
R(1) on Card 10 
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Card Fortran 


JNOe Sis, giNBMes, ete: Description ta 64: Eonmat 

i 13.1 25-36 PHIW(1)) Leave rest of card blank; program E12.4 
j 37-48 PW(1) will determine these inner boundary E12.4 
49-60 SW(1) propertics E12.4 
13.2 1-12 XW(2) Axial distance greater than XLMAX E1l2.4 


(Card 3, Cols. 23-34)(cm). Since the 
inner boundary is an axis the program 
will interpolate linearly between XW(1) 
and XW(2) to get axial locations of 
orthogonal surfaces. Therefore, set 
XW(2) =1 x 10!° 


13-24 RW(2) Final radial position of inner boundary E1z.4 
(cm); since inner boundary is an axis 
set RW(2) equal to R(1) on Card 11 


| 25-36 PHIW(2) Leave rest of card blank; program E1l2.4 

, 37-48  PW(2) will determine these inner boundary El2.4 . 
| 49-60 SW(2) properties EF12.4 

{ 


NOTE: CARD GROUP 14 DEFINES THE POSITION OF THE 
OUTER BOUNDARY (NOZZLE WALL CONTOUR) 


14.1 1-12 XB(1) Initial axial position of nozzle wall (em); E12.4 


equal to X(KMAX) on Card 12.(KMAX-1).1 
13-24 RB(1l) Initial radial position of nozzle wall E12.4 


(cm); equal to R, KMAX) on 
Card 12.(KMAX-1).1 


14,2 1-12 XB(2) Axial position (cm) of second point along El2.4 
| a > a of nozzle wall : 
_. 13-24 RB(2) Radial position (cm) of second point E12.4 
° along nozzle wall 


14,NMAX 1-12 XB(NMAX) Final axial position (ecm) of nozzle wall El2.4 
13-24 RB(NMAX) Final radial position (cm) of nozzle wall E12.4 


{ 

{ 
! NOTE: CARD GROUPS 15-19 AND 20-24 (OR 25-32) ARE 
? INCLUDED ONLY IF PARTICLES ARE IN FLOW, 
: i.e. IPART = 1 (CARD 4, COLS. 46-50) 

| 


15 1-10 FFF Factor which multiplies particle/pas E10.3 
drag coefficient (see Section IC. 3) 
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Card Fortran 
No. Cols. Name Description Format 
15. 11-20 FFG Factor which multiplies particle/gas E10.3 
heat transfer coefficient (see 
Section Il. C.3) 
21-30 CL Liquid particle specific heat (cal/g-°K) E10.3 
31-40 CS Solid particle specific heat (cal/g-°K) 510.3 
41-50 HTRAN Particle heat of solidification (cal/g) E10.3 , 
51-60 WT Particle molecular weight E10.3 
16 1-10 RHSS Particle density (g/cm?) E10.3 
11-20 Not used; leave blank 
21-30 Not used; leave blank 
31-40 TPS Particle solidification temperature (°K) £10,3 
17 1- 5 NPG Number of particle groups (8 maximum) 15 
6-10 NC Index noting whether particle properties I5 


are constant or variable along initial 
orthogonal surface 

0 - constant 

1 - variable 


NOTE: NBL( ) INDICATES RADIAL EXTENT OF EACH 
PARTICLE GROUP ALONG INITIAL ORTHOGONAL 
SURFACE (SEE FIG. 2); DEFINED AS LAST 
STREAMTUBE NUMBER CONTAINING PARTICLES 
PLUS 2. FOR EXAMPLE, IF THE FIRST PARTICLE 
GROUP EXTENDS TO STREAMTUBE NO. 12, 

NBL(1) = 14 


18 1- 5 NBL(1) Last streamtube containing first particle 15 
group along initial orthogonal surface 
(streamtube number + 2) 


NBL(NPG) Last streamtube containing NPG particle 15 
group (streamtube number + 2) 


19 1-10 RP(1) Radius of first particle group (cm) 8E10.3 


RP(NPG) Radius of NPG particle group (cm) 


Card Fortran 
No. Cols. Name Description Format 
NOTE: IF PARTICLE PROPERTIES ALONG INITIAL ORTHOG- 


ONAL SURFACE ARE CONSTANT (NC = 0) CARDS 20- 
24 ARE INPUT; IF NC = 1 THE NEXT CARD GROUP IS 


NO. 25 
20 1-10 WI(1) Velocity of first particle group in stream- E10. 3 
line direction (cm/sec) 
11-20 WI(2) Velocity of second particle group in E10.3 
° ‘ streamline direction (cm/sec) ° 

e * ° e e 

: WI(NPG) Velocity of NPG particle group in stream- E10.3 
line direction (cm/sec) | 
q 21 1-10 = VI(1) Velocity of first particle group normal £E10,3 j 


to streamline (cm/sec) 


11-20  ~=V1(2) Velocity of second particle group normal E10.3 
° ° to streamline (cm/sec) ‘ 
e e e e 


VI(NPG) Velocity of NPG particle group normal to E10,3 
streamline (cm/sec) 


f 
22 1-10 TPI(1) Temperature of first particle group (°K) E10.3 : 
11-20 TPI(2) Temperature of second particle group E10.3 
° e ( K) e e 
° "e e ° 


TPI(NPG) Temperature of NPG particle group (°K) E10.3 


23 1-1° RHP(1) Particle cloud density of second particle E10.3 
group (g/cm’) 
11-20 RHPI(2) Particle cloud density of second particle E10.3 
ia : group (g/cm’) . 


RHPI(NPG) Particle cloud density of NPG particle E10,3 
group (g/cm’) 
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Card 


No. 


24 


26.1 


Cols. 


NOTE: 


NOTE: 


Fortran 
Name Description Format 


IF A PARTICLE GROUP IS AT THE SOLIDIFICATION 
TEMPERATURE THE AMOUNT OF HEAT WHICH HAS 
BEEN TRANSFERRED FROM THE LIQUID PARTICLES 
(i.e. SOME FRACTION OF THE TOTAL HEAT OF 
SOLIDIFICATION) IS INPUT ON CARD 24. IF THE 
PARTICLE TEMPERATURE IS ABOVE OR BELOW THE 
SOLIDIFICATION TEMPERATURE, DENGI(J) = 0 

(J = 1, NPG) 


DENGI(1) Heat transferred from first particle E10.3 
group at solidification temperature (cal) 


DENGI(2) Heat transferred from second particle E10.3 


group at solidification temperature (cal) 
e 


DENGI(NPG) Heat transferred from NPG particle E10.3 
group at solidification temperature (cal) 


IF PARTICLE PROPERTIES ALONG INITIAL ORTHOG- 
ONAL SURFACE ARE VARIABLE (NC = 1) CARDS 25-32 
ARE INPUT (SEE FIG. 2) 


w(1,1) Streamwise velocity of first particle E10.3 
group at r = 0 (cm/sec) 
W(1, 2) Streamwise velocity of second particle E10.3 
° group at r = 0 (cm/sec) ° 
e e e 
W(1,NPG) Streamwise velocity of NPG particle E10.3 
group at r = 0 (cm/sec) 
V(1,1) Normal velocity of first particle group E10.3 
at r = 0 (cm/sec); generally = 0 
V(1, 2) Normal velocity of second particle group E10, 3 
° at r = 0 (cm/sec); generally = 0 ° 
e s e 


V(I.NPG) Normal velocity of NPG particle group E10.3 
at r = 0 (cm/sec); generally = 0 


13 
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Card Fortran 
No, Cols. Name Description Format 
27.1 1-10 TP(1,1) Temperature of first particle group at 10.3 


r = 0 (°K) 


11-20 TP(1, 2) Temperature of second particle group E10. 3 
: at r = 0 (°K) 


1 
e e 


onan, a 
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TP(1,NPG) Temperature of NPG particle group at 10,3 
r = 0 (K) { 


28.1 1-10 RHP(1,1) Particle cloud density of first particle E10.3 
group at r = 0 (g/cm!) 


11-20 RHP(1,2) Particle cloud density of second particle E10.3 
group at r = 0 (g/cm’) 


RHP(1,NPG) Particle cloud density of NPG particle E10,3 
group at r = 0 (g/cm?) 
! 
29.1 1-5 ICOND(1,1) Index which indicates whether first parti- 15 
cle group is at solidification temperature 
atr=0 
ICOND = 0; No 
ails; Yes 


6-10 ICOND(1,2) Same as above for second particle group 15 
e e ® 
e e e 


ICOND(1,NPG) Same as above for NPG particle group 15 


30.1 1-10 DENG(1,1) Heat transferred from first particle E10.3 
group at r = 0 at solidification tempera- 
ture (cal) (See NOTE on page ) 


11-20 DENG(1,2) Same as above for second particle group E10.3 


DENG(1,NPG) Same as above for NPG particle group E10.3 
25.2 
26.2 


27.2 Same as 25.1, 26.1, 27.1, 28.1, 29.1, 


28.2 and 30.1, except for first streamtube 
29.2 ° 
30.2 
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Card Fortran 
No. Cols. Name Description Format 


25.(NBL-1) 
26. (NBL-1) 7 | 
i ae Same as above except for last stream- 
29. (NBL-1) tube containing pa -ticles | 
30.(NBL-1) 


NOTE: CARDS 31 AND 32 LOCATE THE INITIAL BOUNDARY 
OF EACH PARTICLE GROUP 


31 1-10 TRBDY(1) Distance from axis, along initial orthog- F10.3 
onal surface, to boundary of first particle 
group (cm) (see Fig, 2) 
11-20 TRBDY(2) Same as above for seccnd particle group E10.3 | 


TRBDY(NPG) Same as above for NPG particle group E10.3 


32 1-10 DNDSP(1) Distance from outer boundary of last E10.3 
streamtube containing particle to boundary 
of first particle group at initial orthogonal 
surface (cm) (see Fig. 2) 


11-20 O'NDSP(2) Same as above for boundary of second E10.3 
e e particle group ° e 


* _ DNDSP(NPG) Same as above for boundary of NPG 
particle group 


NOTE: THE FOLLOWING CARDS CONTAIN THE REACTION 
MECHANISM AND RATE COEFFICIENTS. USE ONLY 
IF IKINE (CARD 5, COLS, 41-45) IS GREATER THAN 
0. (SEE SECTION I.C.2 FOR REACTICN AND RATE 
COEFFICIENT TYPES ) 


33.1 l- 4 IZD(1) Species A A4 
7 + sign 
8-11 IZD(2) Species B A4 
14 + sign 
15-20 Blank or M A4 
$5 


Card Fortran 
No. Cols. Name Description Format a 
33,1 2i = sign ag 
22-25 IZD(3) —- Species C A4 ; 
28 + sign (if needed) 
29-32 1ZD(4) Species D A4 
35 + sign (if needed) | 
36-39 IZD(5) Species E A4 
49-50 IRR Reaction type (1 to 10) 12 ; 
51 IRT Rate coefficient type (1 to 8) Il 
| 52-59 RC(1) Pre-exponential factor, A E8.2 
(cm-molecule-sec units) 
| 60-63 RC(2) Temperature exponent, N F4.1 
| P 64-72 RC(3) Activation energy, B (cal/mole) F9.1 


33. IKINE Same as above for IKINE reaction 


IV. PRELIMINARY RESULTS 


This section gives the results of several sample calculations made with 
FULLNOZ and an analysis performed to determine heterogeneous electron- 
ion recombination rates in solid propellant nozzle flows. 
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Figure 3. Comparison between Method of Characteristics (MOC) 


and FULLNOZ calculations of MM-Stage 2 nozzle 
pressure distributions, 
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Figure 4. Comparison between Method of Characteristics (MOC) 
and FULLNOZ calculations of MM-Stage 2 nozzle 
gas ternperature distributions. 


A. Sample Calculations 


Initial calculations with FULLNOZ were compared with calculations for 
the Minuteman, Stage 2 nozzle using a two-phase constant y method of 
characteristics (MOC) code.’ The purpose of this comparison was to check 
the numerical accuracy of the code. Input data for the calculations are given 
in Ref. 23. 


Figures 3 and 4 show the pressure and temperature distributions along 
the centerline and wall (with no boundary layer effects), The pressure dis- 
tributions compare very well, but FULLNOZ temperatures are slightly 
higher than those calculated via the MOC code. The centerline gas tempera- 
ture calculated with FULLNOZ shows a more pronounced effect of gas/parti- 


23. Pergament, H.S. and Mikatarian, R.R., "Predictions of Minuteman 
Exhaust Plume Electrical Properties," AeroChen TP-281, July 1972. 
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cle interactions than was demonstrated by the MOC code. The results of a 
short run with FULLNOZ, including a set of 10 reactions involving H,/CO 

chemistry, (see Ref. 23 for the reaction mechanism and rate coefficients) 

are also showr on Fig. 4. 


Figures 5 and 6 show the influence of boundary layer heat transfer and 
shear stress on the temperature and velocity within the wall streamtube. 
The largest effects are observed for the cold (500°K) wall; when the wall is 
near the adiabatic wall temperature .(= 3000°K) the results are similar to 
those for no boundary layer. 
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Figure 5. Influence of nozzle boundary layer on wall streamtube 
tempcrature, 
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Figure 6. Influence of nozzle boundary layer on wall streamtube 
velocity. 


Figure 7 shows the boundary layer velocity and temperature profile at 
the nozzle exit plane; Table | demonstrates that the boundary layer displace- 
ment and momentum thicknesses are much less than the wall streamtube 
thickness, i.e. all boundary layer effects are confined to the wall streamtube. 


Figures 8 and 9 show the resul<s of a parametric series of calculations 
in which the gas/particle drag and heat transfer coefficients were arbitrarily 
varied (via the factors FFF and FFG on Card 15) over their approximate 
ranges of uncertainty to test the effect on exit plane znd particle properties. 
Figure 8 shows that varying FFF and FFG can have significant effects on gas 
ternperatures and velocities. Figure 9 demonstrates that the 4, diameter 
particles can be eitl.cr at the solidification temperature or completely solidified, 
depending on the value chosen for FFG. 
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Figure 7. Boundary layer profiles at MM-Stage 2 nozzle exit plane. 


TABLE 1 


NOZZLE EXIT PLANE BOUNDARY LAYER PARAMETERS 
MM-Stage 2 (Ref. 23) 


Nozzie Exit Radius = 60.9 cm 
Axial Distance From Throat = 134 cm 


Wall Temperature, °K 


500 3000 
Wall Streamtube Width, cm k. 24 1,38 
Boundary Layer Thickness, cm 4,45 2.50 
Displacement Thickness, cm 0.451 0.196 
Momentum Thickness, ¢m 0,355 0.161 
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Figure 9. Influence of particle drag and heat transfer coefficients 
on exit plane particle temperatures. 
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B. Heterogeneous Electron Ion Recombination 


One of the unanswered questions in determining nozzle exit plane electron 
mole fractions is the extent to which heterogeneous electron/ion recombination 
on the surface of solid particles can enhance homogeneous electron/ion recom- 
bination rates. Towards this end we have adopted a model developed at 
AeroChem by Calcote, Kurzius and Silla*‘ in which a negatively charged parti- 
cle (particles are negatively charged because the electrons will reach the sur- 
face more rapidly than positive ions) is neutralized by positive ions striking 
the surface. Thus. instead of requiring a three-body collision for recombina- 
tion, only a two-body collision between the solid particle and positive ion need 
occur. Toa first approximation the particle electron recombination rate can 
be equated to the rate at which positive ions strike the solid particles. How- 
ever, if negative ions are present (for typical solid propellants mole fractions 
of Cl” are from 2 to 3 orders of magnitude greater than electron mole frac- 
tions) the electron recombination rate will decrease, since some of the posi- 
tive ions will react with negative ions rather than electrons. 


The solid particle electron recombination coefficient, Ope is defined from 


dn 
© ae 
(GE), =-ape BH Mp (71) 


where ne is electron density, t is time, ny is positive ion density and Np is the 


particle mumbo s density. If particle diameters are small compared to the gas 
mean free path’ (free molecular flow) and the electron densities are suffi- 
ciently high that the particles remain negatively charged, ap. is essentially 


the random ion flux to the particle (with a correction factor for negative ions), 


-1 


Y, n m_ | 
ma) ‘(Se a (72) 


t 


If particle diameters are not small compared to the mean free path correc- 
tions will have to be made to the expression fora. 

e 
24. Calcote, H.F., Kurzius, S.C., and Silla, H., "Solid Propellant Flame 
Ionization and the Effect of Chemical Additives,'' Third Radar Attenua- 
tion Symposium, CPIA Fubl. No. 46 (Applied Physics Lab., Johns 
Hopkins Univ., Silver Spring, 1964), pp. 17-40. 
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where r is the particle radius, k is the Boltzmann constant, T is the gas 
temperature, m,,mMe, and m_. are the mnasses of positive ions, electrons and 
negative ions, respectively, aid n_ is the negative ion density. 


Equations (71) and (72) represent the formal method for incorporating 
heterogencous electron/ion recombination into FULLNOZ. However, we must 
also account for the possibility that, at the high temperatures near the nozzle 
throat, the particle is emitting electrons via thermionic emission. Under 
these conditions a steady state is achieved by balancing the random current 
density to the particle (from the plasma) by thermionic emission. Equating 
these electron currents results in a "critical" electron density for which the 
net current flow to the particle is zero and an initially neutral particle will 
remain neutral. This critical electron density is defined by, 


(nglee = B (#) 7, ex| ~11, 605 Eq/Tp| (73) 


where B is the thermionic emission constant, T, is the particle temperature, 
T, is the gas temperature and E,, is the effective work function of the particle 
in volts. Thus, heterogeneous electron/ion recombination is only of potential 
importance for neg > (ne)oy- 


The technique adopted to incorporate the above equations into FULLNOZ 
is: 


1. At each integration step determine whether the local electron 
density is greater than the critical electron density. (During the 
initial stages of the expansion process, where particle tempera- 
tures are very high it is likely that ng < (ng)cy-) 


2. If ne < (ng), then electron/ion recombination will not be significant 
and Eqs. (71) and (72) will not be employed. 


3% If ng > (ng), Eqs. (71) and (72) will be incorporated directly into 
the general kinetic scheme (with possible corrections to Eq. (72) 


due to non-free molecular flow effects), 


The above procedure has not as yet been incorporated into the code. 


55 


10. 


ll. 


V. REFERENCES 


"ICRPG Two-Dimensional Kinetic (TDK) Nozzle Analysis Computer Pro- 
gram," Dynamic Science Corp., December 1973 (revised version). 


Hoffman, R.J., English, W.D., Ocding, R.G., and Webber, W.T., 
"Plume Contamination Effects Prediction: The CONTAM Computer Pro- 
gram," Final Report, Air Force Rocket Propulsion Laboratory, AFRPL- 
TR-71-109, December 1971. 


. Nickerson, G.R. and Kliegel, J.R., "Axisymmetric Two-Phase Perfect 


Gas Performance Program," TRW Systems Report No. 02874-6006-R000, 
Vols. land I, April 1967. 


Penn, M.M. and Smith, S.D., "Supersonic Gas-Particle Flows, Includ- 
ing Reacting Chemistry, '' JANNAF 8th Plume Technology Meeting, 
Colorado Springs, CO, July 1974. 


Boynton, F.P., 'The MULTITUBE Supersonic Flow Computer Code, "' 
General Dynamics/Convair GDC-DBB 67-003, February 1967. 


Boynton, F.P. and Thomson, A., "Numerical Computation of Steady, 
Supersonic, Two-Dimensional Gas Flow in Natural Coordinates, saa i 
Computational Phys. 3, 379-398 (1969). 


Van Driest, E.R., "Turbulent Boundary Layer in Compressible Fluids," 
J. Aeron, Sci. 18, 145-160 (1951). 


Keener, E.R. and Hopkins, E.J., "Van Driest Generalization Applied to 
Turbulent Skin Friction and Velocity Profiles Measured on the Wall ofa 
Mach 7.4 Wind Tunnel,"' AIAA J, 11, 1784-1785 (1973). 


. Pergament, H.S. and Kelly, J.T., "A Fully-Coupled Underexpanded 


Rocket Plume Program. (The AIPP Code). Part I. Analytical and 
Numerical Techniques, '' AeroChem TP-302a, AFPRL-TR-74-59, 
November 1974, 


Pergament, H.S. and Kelly, J.T., "A Fully-Coupled Underexpanded 
Rocket Plume Program (The AIPP Code). Part UI. Program User's 
Manual, "AeroChemTP- (in preparation). 


Mikatarian, R.R., Kau, C.J., and Pergament, H.S., "A Fast Computer 
Prograin for Nonequilibrium Rocket Plume Predictions, '' Final Report, 
AeroChem TP-282, AFRPL-TR-72-94, August 1972. 


12. 


13. 


14, 


15. 


16. 


17. 


18. 


19, 


20. 


ra 


22. 


23. 


24. 


Cruise, D.R., "Information Manual for the Theorctical Propellant 
Evaluation Program, '' Naval Weapons Center PEP NOTE TN-U-1 (plus 
additions), December 1964. 


JANNAF Thermochemical Tables (Dow Chemical Company, Midland, 
Mich. ), continuously updated. 


Jensen, D.E. and Jones, G.A., "Gas-Phase Reaction Rate Coefficients 


for Rocketry Applications, '' Rocket Propulsion Establishment Technical 
‘Report No. 71/9, October 1971. 


Soo, S.L., Fluid Dynamics of Multiphase Systems (Blaisdell Publ. Co, 
Waltham, Mass., 1967). 


Crowe, C.T., ''On the Momentum and Heat Transfer Equations for Two- 
Phase Plumes,'' Washington State Univ., March 1971. 


Schlicting, H., Boundary Layer Theory, 6th Ed. (McGraw-Hill, New 
York, 1968), p. 599. 


Herring, H.J. and Mellor, G.L., "A Method of Calculating Compressible 
Turbulent Boundary Layers,'' NASA CR-1144, September 1968. 


Tong, H., Buckingham, A.C., and Morse, H.L., "Nonequilibrium 
Chemistry Boundary Layer Integral Matrix Procedure," Aerotherm 
Final Report No. 73-67, July 1973. 


Hopkins, E.J., Keener, E.R., and Louie, P.T., "Direct Measurements 
of Turbulent Skin Friction on a Nonadiabatic Flat Plate at Mach Number 
6.5 and Comparisons with Eight Theories,'' NASA TN D-5675, February 
1970. 


Back, L.H. and Cuffel, R.F., "Relationship Between Temperature and 
Velocity Profiles in a Turbulent Boundary Layer along a Supersonic 
Nozzle with Heat Transfer," ALAA J. 8, 2066-2069 (1970). 


Kays, W.M., Convective Heat and Mass Transfer (McGraw-Hill, New 
York, 1966). 


Pergament, H.S. and Mikatarian, R.R., "Predictions of Minuteman 
Exhaust Plume Electrical Properties,'' AeroChem TP-281, July 1972. 


Calcote, H.F., Kurzius, S.C., and Silla, H., "Solid Propellant Flame 
Ionization and the Effect of Chemical Additives,!' Third Radar Attenua- 
tion Symposium, CPIA Publ. No. 46 (Adplicd Fhysics Lai. , Johns 
Hopkins Umv., Silver Spring, 1904), pp. 17-40. 


57 


APPENDIX A 


SAMPLE INPUT DATA : 


ac a alia ts 


Oct Glen ‘ FUMTaam toe foe 


FULL NOE , 
CARD 
v45) - 
s- ! : : 4 , 
2 wes FULLY COYPLLO NUZZLE Slow PROG#AM — fULLHOZ - SAMPLE TEST CASE 9/75 ¥9H 
3 7200 2. 000E+2 1.524842 
4 / / (e) / 0 
5 15 b 10 2 22 40 $0 900 10 
¢ /000° ! of ; 
aa 0°68, 0.0! 1°OL-4 1-0 O-f 
811+ 64608562 E+1 + 1P$75287 Erl-*F¢4IBCSBSELIO' 3628095) £-! co! 
2-590 4g64 £-2-' ABddSSOIE+2Z' 544./0327. Er2 CO 2 
821 sa7g/ SSIES 1 I SS92ZFOGE +02 = Bb056609LI *LISOOAITEFS cot / 
2- $056 7069L 11 - GOFIIOBIZE FL *ShIQDIGGIL GE cone 
B31 49685629 £t1-*23301L6BE-3 S96 67528. £-4-' SB V2 545£°5 wl 
LH /B4467323LE-4° 506/86 S55 Er2 33404654 Er2 H2 
8.4/-60596451 Ler I34/2152 Ee l-+16S5086752REtO*BB760457 E-2 Wit 
hi hhc batt Er /-* /T00SOLIE TL: SB4F43527/. Erz HZ 2 
85) +58018265 Eth 17395975 Et) ~*325362/B£10°22313547 E-/ ‘On ¢ 
2751795 $) £-1°79212009 Ets 50930996 Lt2 on & 
861 +6031 1014 Er)» 47049758 Et) ~*9b629) L45E 10+ 49087996 £-I #20 I 
2G OdZ/IS2 £-1~*S9SD3S34Lt2°SIZBTLES Lz eRe oe ne 420.2. 
LULL POR 2 /# 
CARD a 
” er . : . 
B.7/ 62693469 Le) 193920092 Erl-+46/58274£+0° 37003135 £-! . W2it 
2° 9/GS5276B8 L£-2-°/924¢240B8E4/°S282477¢ Led Wz 2 
BBI 50617760 Etl- "17366 SEEtO'394053$8 £-/-*2O83IS77E-2 Oo! 
(2° 18127099 E-] 58115304 Led 447L8/72 L£t2 oO 2 
B26 +724/7509 br) (2398/22 Erl-*1B448469£+0° 10507580 £-/ on 
2-596 0406TE-1-' 24219 96BE+1*57031740 Et2 02-2 
Bia! -4B4B5464 Er] 32006164 £e0-* 2054547 /£'0'BS218696 £-I kl 
(2°43725018 £°2'19867370 Leh 441065/4 Ltd K2 
9' CO /39+ E-6 300° O75 /:¢ If 28° 
2602 {39° £°6 300° 0°75 [4 1-4 eee 
3H 995 £°6 500° 0:75 /-¢ 4 fe 
4H2 63'5 Eb 300: 0:75 14 14 2: 
50H [25°85 E-k 300° 0:75 4 /-4 17° 
6H2O $25°5 E-6 300° 0:75 /:€4 4 18° 
TNL /60° £-6 300 0°75 leg 14 28° 
BO (66> £-6 | _ 300: O75 /'4 If /6° 
9 O2- 189° Eo | 300° 0:75 14 P¢ 32° 
10K {29°77 E-6 300° . 0°75 4 /-¢ 39°) 
~NOTE~ Cao 10 Omi7TTED FoR ITURB*O(CARD 4, COL. 56-60) 


Ob Kn ee FORTRAN Colng Fore ; ; 
777) Wee ae ane foal 
Re Se te eae oe rete To ge A Ngee a as 
i” Z aa, - iar, a 7 7 “Woetaans startet ioe . 
See 2s eee saa ee 
i eae "ole ovao¥  gls00 
hae ‘BE-4: 2 #£-2 2, 66-3: i, oe 
he | Hy H 
htt | Te. 1h. bt:2| ALE VOENTYCAL To ©: 
: 22 eg rs 7 
= oe ht 
; we Pid PpReesa 
\ i \ | i | L ie ee 
CARO ; ; — FO TAN, HATESIAT S is eras a 
271 Ve Aol. i deere | 6, eae 7 | : ig ee ey ar Sshed 
, i | | 
ie ese ats We Gb a CoP wa | 
128! ‘12-319 48278 0-14428 | B+ 103 2739", 1+ S50EtS 
yi (SER OTE a ! 
a is ‘ | | | 
29). a-aat | 5 ROO) 16232 © | B+ 003 | 2732" I FELE+S| 
3 SER MIE | 
1 | 
ol, ,f2 3; boas: 0,-1.8 0315 ete 2726"! /\ 5626 +5, 
1 Se NOTE ‘ ’ ' | ’ ’ | 
: tes : ’ ‘ ' 4 ' 1 ‘ ' 
ani '}a-ops |b b176 Ces 903i9 | tas | 2719. 1 SI0ES 
; 1 i i 
‘ ., SE ed ' { ! ' { 
| i | ip pe | 
t i } ' ' 
eit! 1/8 5 | 72100 | ‘ 2164 | 743 |! e7/2°) /i $76L+5 
é 1+ SER NOTE ; . 
3, ’ ! : 4 ' | ’ r ' ‘ 
! ! i Lich &. o.0 Z 5 a a4 


FORTRAN Loading Fura 


CAD 
No 


vei ‘Vora i gero.” “O-S287 Fee ges 2594+ fb TGESS 
¢! 5 SE! WOTR | By dl 


aol "10° 507 11+B500. 036070 5+ 950 25e3° 
' SEE: NOTE ; ! 


"TLSE*S 


' ' i ' 


zs ~~ 7 RAR Aelita Fe __.~nneeeneenrei 
4 Ob Fin 2 FORTRAN Coding For : ; 
. . . FULNOE Ont eee ores) On, a bllos.  Bie ee ee 
4 this sey oak 5 = ae ae oe a 
; hi . _ a 3 ie , iz - ; . -* 7 x S stgamaa “tastes a f ? ~ 
: [7-000 | O-23480 7-98. [2708-1 S88e+5 - 
3 a aoe 6 ae ns & nd | 
4 2 ivy 2524) | gel 2700+ 1 600E+5, 
‘ 2 ' : i i \ 
_ ee ema ob cme x ats or P H } ; 1 aes ‘ 
sh pee eg Shee teal aes ae ee 
3 : 1. &+|97. 70513 “9 | 26 88 «| (6 19L 45 , <a 
4 Zi i ieget, | aaa ! 
ce | 
4 ss EZ 9-546 |” een ew 0 Ye ee 267/+| Ire aero. | 
q i i 4 
3 : ; 3 Sdeeee es ' - { 
a aay box i ' ‘ t { ‘ 
: , i 7 I ‘ | 2648", 1647 b£t5, | 
~. {Worsia 7-090 | oe saeels | BO | peed] fk 70E+F 
a ; 2 ete ‘ ' ! ' , : 
7 bee, Ua : aes a -- | ‘ + . t 4 
: - ete iat ee teh ect ok a - at oe te RS, Mea parc eas aiken ae ~ nnn 
: 22 pone | a641$0° 037874) | Seka 8550 I REFS. 
d os SEE VOTE 
2a). 10° O60 /2°9750, 039677 | 5390 2533° «| 79BEFS 
i a SEE: MOTE : eo Ns . 


4 

PI . 

t RBs 9 BLDS peasy O 4148) . $065 2562° 1 847895 
s 


ee 1: StE;AOIG . ‘ ; , : 
He t ‘ . . t ' { 


9-575 14 0660 Faces oz | 2401+ 1° 904E+5 
LE WOTE 


OB 6-Pen : ae FORTRAN Codey Ferm : 
Do Fueenoz | S , ge to I : ; ' 7 ¢# | 
me i 1 H i d 
— oO) aa oe a ro Y eat T 1 i 
i 0:0 I = | | | | : ' t 
/4:0B84 : 
: 140275 | { 
| Mesos; jy Eg | 
' L521 | i : . ; : 
ay el i | | | 
: (6:°434 Poe ; : 
17149 | | i 
i 18999 ; ' | j { 
202/68 | | 
22-606 j 
24|-994 | | 4 | | 
27| 305 | ; i 
30\'683 | | | | 4 
331-985 { : ; ; j 
37; 084 | i 
42/062 re a ee ae 
_ 44838 t | | | | j 
50) 149 | | ' 
54/508 ee | 
$8039 ; 
6 0\909 | | ‘ 
ESO ns nee eee a een ee | 
FORTRAN Coding Form 
: Pe RO 
! ir oe 
Vd" 3st ~ prse assy Loom | 
BO, ne he i 
| | : : : | i | 
i are j 
o t-OE-4 3:0£-4, 4 OE-4 ' 
ly WC peo 17, cof. 10)! 1, THEREFORE CARD NO 25 1S. WEXT ; 
LS! 0°14 Erb 0-139 Leb. 0-13) Erb, 0-124 Ene . 
261 -D' 27h, E¢3,0° 559 £+3-0° 768 Er3,-0-/0) Ee 


2EI B13 2893 29/3° 2973+! 

2B 0 279, E40 143 E-3,0°/82 £3, 0-554 L-F 

£91, OO, 0, O} ; | 

B01) | 00, 0:0 i: iD. | 0-0 | 
‘ i 


25.2 0 IAT Erb, 0+ 139 heb, 01/32 Erb 0-125 Eve 
Ui-9- We £03, 0: S92 Le3-0°/4 Ed, 0-107 Led 
oe £8/ 3°, 2899" 2x0 ds 2975" ; ; 
parr é # 0° i$ E-3 0+ /82 E-3| 0- 5516 E- 
ry, CARDS 29.2 T029.18 fee 10k "CA 70 29.1 ' 
| IO! fm ONL CARDS 30.2 TO 30:/8 ARS HOENTICAK TO 301 ; 


4 
263 p-s4B Soil o: / 39 kee 0-/32 £46. O+ 125 £6) ; 
263 D'3/p, EF 3i-0°62P £43-0- Beg £t3-0- im Eth 


ws 40-1 OS: 

| Et 30° 986 £+ 
fara alae 
28. ee Ord: . E 


iz ft a 


FORTRAN Coding Form 


~~ Lk MOE 

Uk Dead cree gee een ite,” ee | 
caeony 

es 


306 


FORTEAN TATEMENT 


ance! 


sae = aces: oe ha 

wert p- 14 pil 140 E+6| O° F evel 0°12)5 E+6 

ue 7-(0 4114, £+3- 0-839 L£t3-0-/15 E+4-0-°151 E+ 

an7 | 28013" 2003, 29138 2963° 

28.7, D279 EM 0+] 43 £3, 0/82 a O° 554 
1 i { ' i 

a a 

43.8 0+ 14,9, Erb, OF) Ff £66. 0-133 Erb. O° val Lr6 

Ree AL Ee O'91K £43-0°/ 24 Et9\- O° 165 Erd 


| 


278. 280)" | 2883: 29/3" 2963") | 
2B. p-278, £ a4 143 £73, U8 £-3, 0-554 EA 
308, I! 
as | ; : , ; 
259, p-/ 5. Evel O° /4) £ +6, 0+ / 33 £+6. oO 126 Ere 
2L.9,-0° 499 Et5|-0+ 10} Erg “0°13 Led: 0 a Er¢ 
£79 | 2B 0)3" 2693)" | 291 3" 963° 

89 pap EA 0-149 £73] 0-182 E- 3, O° 55 £4 
£29 | i ; ! 

30 | . 


| ; 
t ’ ; ’ . 
he | } H s ! 


Dob Fn FORTAAN Coding Porm : 
BE AMEEN ORs oh Gen, caeth eat any aaa ae aes a Ms (oa 2 &, SR 
' | } [ a ee Be ee wee a tates 
By ST) 
ees Wheat) Juteees amicnane: eee eee 
! i ‘ 
| ; | | 
[. ee 
i paTm c-4 F149 Ea 0° : 
aah) | : 


20 a ae 


si) F 4 Ps Ov! ‘ 


26-1 + 624! £43| O° 73. 


: 
a he 
’ = ne EON ee eee 
Bit cousiee wnelsSeae Mates anni les Eats: EES 


3 
2 
z 
s 
3 
opens tvantcatle ictent ta Paihia us atin Rett wad Enn ven tniaart Ste Sebo icktsle eS 


° Pe £+6) 9 rf Fes 
73 E£t#h-' zr Er: 


, 4 : 
Se els ORIN REO 


t 
: | 
iy 
281", £4 0+/48 £-3| 0° S¢ - | 
i 3 3 ae s 75 4 Eg 1 { 4 
2917 be? 93 ‘ | : 4 
is eee te be ie ) en ee aa ee ee 
‘oe Roto 4 oe % | ‘ 

2.2) 3 Ercl ov gle Evel O-saiy Erc, | | 
4M: bony Erf-0- 196 £t#\-b°260 Er ; ; i 
2 | 2760! 2873)" 2893) 943° | | 
$l l0-n719, E- 40°49 E-3| O° 18R E- 31 DSS EA | ~ = “abe 3 : 
A322 es eae, sree” 1 sf 
a ee ae 
25/3) sl god Ol bt Ere) O- 13h Ets a0: 126 Ete | | | 
uss a Soe | 


FULNOE ; ts a Ge che ah pten oh ee RAE OIES kata 

i ‘ bok <<: ceiy 

Oe sii ' a 

Be oa78 E4 0-143 E35 0- joa ‘E- 3] 0-554 é-¢ oe a cas ake ates a a 
j H : ‘ ' \ i ‘ \ | ‘ 
ae ae ae | Ye Yl ! | 


U4 0-154. Ebb Of 4S Eth O°/ 36 Ete O° /1B Ere, 
weit-o- 993 E+ 3-020! Et#-O1277 EFb-O-363 EFG, 


4 2763 | 2883) pees 2939° 

BM 279. E-4 O14 E- 30 JOR 2-3 0-554 E-¢ . 

he ae ae be Fai si ie ge | | 
30/48 } : 
BIS OFF E+6l 0-146 ie 0°/37 Ere, 0/29 P| i 

US -0° 12S Ergi-O-252 Et h-0-34 Erd-o-45t E+4 ; | 

2718: 2763 "a 883° 93 3°! H : 

WIS p:279 E- O+ 143 E- 3) 0+/82 E-3| 0554 E-4 | 

29S _ ' : ; : : ; 

ae ; ; ; 

£54) 0+) Sp: £ +6) 0-/47 E+ O° / 38 ay ote Eee | | f : | ; 

Loe 10+1 ob: Erd-0°3316 Erd-Ordbi Erh-O°boe Er4 | : ; 

ah a7 4 [ 2623) arse, | 29239) | : 

hy: is E- ‘ 0-143 E-3,0°/ 02 £-3, 0° 554 E-4 

29% ; ; ; ! i ; \ 

nc ne ore ; i | : 


ObeFhen FORTRAN Loding form 
777.7) Se | 


a S Beast $ >). ae 
Dit, ete 


tee TATE MERE 


Se ee ea 


Pee pT a. 


oe craton 2itiete Sen emiieeee. - 


neo os Ee Sa ETE EE see oe 


BI ER Pees pore 4 il kh Ot 


i 
i 


er a et ere cin 


APPENDIX B 
SAMPLE OUTPUT 


COLEY COUPLER ZZLe FLU® HRUGREO 


on, a 
ceaae HER AN 


nse SiR stdin seed 


AERUCHEM RESEANCH LAGUWATONIES, INC. 
PUINCETON, ake JERSEY 00500 


t 00 FULLT LINPLED MUZZLE FLUP PnuGHA © FULLAUZ © SAMPLE TEST CASE O775 ene , 
a tte to date Sey 1 
sat HSES Fevu VALDSS Punk ove atebas 1.9200 002 aL“yosee, Imus oe ants = 38 $ 
{ 
3 4 
4 Aman fn a w ict n\n )) Ty) i) Wile 10106 = LPL ant i 
3 o 6 te C) 1 ' é a 1¢ me se e 906 : 
aw tele PCr tome leuEse = Fte 1csn BraPes inine = iPANt ing itues = 1P te 4 
3 oo «. en Q ey 6 ) 1e ' of ry 6 4 
A : 
% ne tae ww mln : 
‘4 t,visut ove Setezut eos 1. PHiteny 3 
3 FA 
3 Auenan — eesuUN fuL OReiA atm tater GAav banc Fuactn y 
t oAue us Dey Cu oe tou otud 00,60 0.00 6.06 $0.00 4 
; i 
ry tote ‘uv 0 UML GA on ac “a ‘ 
FA ; cu te s¥uteys 5.000t 02 Ted80G 001 2.800E 000 taeoteoe i 
q Cue 1, S¥utede 3,000 902 7eS@0E $04 $.000b000 i 
8 " 4, Souteud $.00ut 002 TedOCE°O1 1.600080 : 
3 We a, Sduf oud SovCueene 7.Suvteus 1.G00E 60 2.0008 000 Q 
¢ un teedrf ote evun eee 7,S00t ous 1eaudeeoo 1,706 061 4 
4 neo teed%h ous S.u0tt eve 7.5u0E ous 1.veeeus 1.800t 001 , 
4 ne teouvéeua $.000¢ 002 7.S00t0! 1.aneceoo 1,aceteoe 2.000000 4 
4 ) Leondceos Se000E eve 7.900803 1.400E cee 1.400800 1.000800) i 
i ue 1.890E°0u 3.0006 ¢02 7,S00tus 1.4006 000 1. eget eee 3.2008 001 
iH « 1.297E 000 3.00UE 002 7.500F on! 1.4006 00 124006000 4.0008 008 3 
‘ io 
4 «oa w lad | v ' au tT] : 
2 t tedoteoul Ue 0, : i 
we 2 tedooteut @,vssteus 1,boateoo2 8,S00Ee00 2.7etoos 1,.53ete0s H 
zal 3 tedesteor Teerte evo 3,0u7€e02 6,87 3E000 2-759€003 10930E 905 
. @ t,Pedteor te Blot ene SeG1tF ove H,aa3Eeoo 2.758E003 $.535£ 005 : 
i % beddttoul Se GehE ony Deelatoud 6.800b 000 27ST ONS TeSS7E COS 
& Leedstout Seveatevu O01 6b one &, 30$€ o0u ee Tht ous 125996 40% i 
? teeuteovt S.odotouy Levh2t oo H.273t 000 2oPdicous s daltons j 
@ teddutous 8 .P27E 400 1l.2edien) H.19SEC09 20 78SE003 1.5406 +05 ‘ 
© W.252teus 4 .A2Ate00 1,443¢00! 8. 108E000 207396003 1.Soseees : 
10 1.22Stony Seueot out 1,e23t eu) Lube a 2.7326 o0§ 1 .550F 405 4 
it LePtetou bece shou bebuat eon} 7, 8956 ¢00 2eP2ot 03 teSesteos } 
12 beeuteour BoOtme sue 1, 9ret ott 7e77Se 000 eo71Geeus 9.5716 405 5 « 
13) 1 tevteut fePiut oun eohoutous Teoetiouu dob teteods L.d79be05 + 4 
te F,t7Steouy T.MOUE Con 2. 30SE 004 7 ,Q@8E 0" 2,7088¢03 1o580E 00S . ‘ 
1S 1 dovesot 6, 388¢ 400 2.525t 003 1, 3208000 2. 7008003 1.0000 965 : 
te t,sadbour HOTS Onn eo POS oM1 7,1 386000 ee OANE OOS teolateds f 
17 Letenboor Vodhoe een 2oHAot oot ©. 95St U0 2ob71E005 1.030005 4 
18 Uedbueous Lovdseout S.tont Ot eT IStouu d,0ust 003 beheutoud ; 
19 L.O%teout Pouriteur Se2doteout 6, 480E 000 eooeuceus toorikeus 4 
20 1,07ee oor Leldorevs 3.4276 03 ho2eot. 400 2 od%GE 003 1.0706 00S q 
21 teodteeut telasteoot S.ou7t 08 59806000 2.S0sb003 1e720€ eos : 
22 beveveent Tecuet ent 3, /MIE oUt 56806000 25906005 1. 702beos A 
23 b,uvobout Leevek ong SoVont ous 9. 30nEoua PeDS3 ONS 1. 7emtegy : 
26 Vhelr edn PeAdteent Gelube et 5 .06SE oun eo Svet ous 1.h@7b 90S H 
2> WSlet oun tune oot G. Senteout d.nudteouy e.eureous 1, 90st cus F 
{xe wa belie vr Sa 4 
1 Vednleour Ve o, oe Ue j 
@ tetuctotu OU. eon, ou, ov, ‘ 
NS Ae cry vrib pa au j 
1% dfetone 1, dee ay 6, 0, 90. £ : 
2 FNS 7Esn0 t,veatoul 20. oU, 206 ' 
3 FeOSdbeny 14508 ont 20. 0, =, 4 
oe Detooreny Lede2rout Ue ot, ou, M 
S beedmorent ooh dt oul a0 Ne “0, 
mh bemddr ont Vemsse ett PVs om, wlte 
7 Lebaot out 1,27 3e ont se o6, 0, 
A 1,9%eutoul Pawo ent “U. Ue 20. 
Y debfureus evlecteul 20. we Ue 
fu eenbetouy Coeoiteul oe one Ue 
Hi betvurery Peueve out en, OU. oo, 
$2 S.p9nreut 2ofsitont “Ue Us "060 
13 @,seuteut SeVORE ont "0, ed, 20, 
16 Seeeet ony 5. 590k OM} °0, *0. eu, 
> SS Vhageut Sedume ety Sie 2Ue Ve 
fe @,s¥beout wer teteus Ue OU, en, 
V2 mrt teeoul enndd ont , ev, eVe 
te ¥,9usreus wevTSeent “0. 6, 20. 
19 ttbeteve Swit out eb, ev, 0, 
OU Leesoe eu S.muee ous ou, 20. 70. 
Ot by tear ene meevteend PUe 6, Ve 


Pe lWeSdut eur Herm de ord 20. ou, 20, 


wa 
q 


ne 


e 78 : 
Raa tien tat ete me Lame abe vielen. 


as 1, goestent 
soe gleticeus 
Cir mee PE a) Pd 


arecit Cu 
aveCie Lud 
eel on 
areche ne 
wMecit in 
SeCit Hdu 
the = na 
secie 
eeeche ve 
Seecte » 


BS tadosetees 
mae 4127 8e Ces 
vis N7heceoe 


aeecit tu 
oreCie bd 
seecie 
seecle m2 
SrECIE unt 
Brel meu 
Seecie ne 
Seecit uv 
MECIE ba 
SVECIE a 


BS Lecub¥Eens 
os oh2hoeeny 
Ole ,b7upeete 


Seecte te 
Sele bie 
Seecit on 
aweele he 
SrECie Un 
BrECIe Han 
Specie Ave 
Sweetie ov 
wea te 
sPrche 8 


as 1,eonie eu) 
ras oferuceny 
PIs .leUneoue 


seer ie bo 
SPECKE Chie 
SPeChe om 
avin me 
HeoChe te 
weit nev 
SPECK ee 
aeeCle 
SPL le ve 
Le 


aa t,asddeee) 
wae ,1de seeds 
wfa 4 tnvarede 


sercie oe 
S¥kEbe bud 
aweCle 
settle me 
SheCle iim 
a YY 
nee lie 2? 
srevde t 
Srecie ve 
srecie © 


aa Neeaseeort 
AS ghetor ent 
via glenucend 


APPLE oe 
Sorvude lee 
shedab ¢ 
Srelat we 
eeecie ta 
SPELEE meu 
SPrcit se 
seette 
aPoC te ota 
omec de ow 


He ©,.030 00) 
UrLTs 


one 


ce 
ts 
Cs 
ts 
ts 
te 
ce 
cs 
cs 
ts 


WB Leet ioteue 


vELYe 
bhes ‘ad 


oOM4 eve ne 
oJuaStces Sas 


Pugne 6. 
Ouse 31,0000 ce? 


oSOnteces 
Ge 


o3eeet cece Ss .28aseone 
000006 00) BS = 8700G08) 
o&$00t 263 BO 820% 0Os 
oP BUG oat Re 42seogean 
efoentoue ae 4S0S3 ove 
03 9008 066 ae .2e03eeus 


oG20Or eve Be .d¥7SECOO 
evetous ae _veoteose 
00F 206 Be, $e7lgoee 


Siieetous ke NaS atowe 


2hOS0E0CO he 
eFazateods Sus 


o S0UNF Cue xe 


Onis 3,007e 02 


o8OS7ECOS 
% 


o212deeee 


oOb0t oo) Re 27 CGROUS 


e 

of POhh o0u as 
42008 000 as 
oT B00tooe as 
eSOnng o06 xe 
o1fvoeeds se 


WE TeALSME COU vais 


verrs 
nds 


Cs 
Ca 
(s 
ts 
Cs 
Ca 
Ce 
ts 
{s 
ts 


RB SQGIVNr etn 
veers 


knos 


ts 
¢s 
Cs 
ts 
(a 
ts 
Cs 
Cs 
ts 
Cs 


e00Ses ene we 
elauitevs Sx 


o SHOE ODO is 
oO Cys as 
emdulteous as 
oRSUUE COS as 
o2o0vgoud as 
otMoteue ae 
e4PN0t COU ¥e 
o THULE oue RB 
ehrnve ety as 
ethucbeous ae 


eOUSTF OLY HO 
oT 3netons Sne 


08280061 
023006000 
oSuSSECC2 
o2u93E eee 
abPISE ove 
oVeodtous 
9 3071kege 
oD450E CRO 


Swi tueoue 
019526 003 
¢. 


21 25E 000 
of 0GF ous 
ote4%eoo) 
o2580C000 
05035002 
02093008 
0297S 000 
oVaenteys 
oSeliteye 
oSddbEode 


HTS T,a1eueeue 


otoareeos 
, 


oSG0CE C00 AB .8tede ooo 
ebCUUE OU am 2708 oys 
ee S0CF oud RB , 126% ous 
eeaiugeons as ,236uEe90 
ee Orns one se, S9SSFOUe 
eh ¥ungevo Re .2095E 008 


8200 600 ae 4aO7SECU0 
o TANK OUS 18 ,VeedEo 
ehOultone ke =, 3a7lte 


eh OGrous az 


ke 8, Ne $heeoe vais 


Otis 
ante 


ts 
Cs 
Cs 
is 
(s 
is 
ts 
(se 
Cs 
us 


eOOSGE ONY Ne 


aTS2ibous 8x8 

eSuutbeun as 
eOCUVECU) is 
oe 300b 003 is 


Pr as 
ePADLE CUP ae 


ot4ott out as 
ee Fhlk out as 
o fBLUb eba as 
edbnuteene Xe 


of 00G dS xs 


We beocede enn onis 


veers 
ants 


(s 
ts 
is 
ts 
is 
Cs 
ts 
te 
Cs 
ts 


ett Saeeen ne 


POL) 


9,0: 00Eeue 
(057K 003 
0. 


Ph 2oFeuy 
edi Oateout 
of 2a%Eont 
2 SF0E 000 
o W083E C02 
e2u9 st ouy 
eeVsStoue 
oPaottous 
dai koue 
oSeoetoue 


1evddtt oot 
ebthingeos 


eleTthous San 0, 


2 SUM UU as 
ehuuut et a8 
et SUF oH§ 8e 
eeduutel) as 
ePOutd o0@ as 


2B 9U0E 000 is 


eG200F O08 ue 
Pee xs 
ode peve as 
etiuuters 1s 


e2t2SFeou 
ePlvabent 
ot eevteut 
of ShUF Oued 
SUSSE S02 
od09Steuu 
o79TSEeOu 
eVoont ous 
otuliteta 
oPOhab ove 


fe 2.Feteteed 
wis =, weet ees 
SUmouts 


a errata tmenmameme: mai tnn ot sine tne Dil, eh a we 


ob Stab 005 


te 2,.75000003 
wie §=,38ritees 


eunoutes 


oS8B96 003 


is 2, TomvEeud 


nls ,38engeus «= Tans 6, 
BMOUIS .s1ekeee 

ALUTE 0, 

AdliTs 0, 

abuls 0. 


evuts 
avuls 


te deTeivecod 
mta =, seevrees 


susuuts 


suUuts 
evuls 
rouls 
evuls 
mids 
auyis 
evuts 
adultes 
avuts 
evils 


fa 2, 7SSvboo3 
nts ,$Rareees 


SumuuTs 


auuts 
euvls 
avots 
aunts 
avuts 
bots 
auits 
avuls 
abuts 
avats 


Ts go fSiutess 
mis ,3A87bCUd 


Susvuis 


avUs 
aes 
aunts 
ails 
wuts 
avute 
CU 
etoile 
aude 
avuls 


80S tees 


Tans 0, 
oSdSoccee 


aes 0, 
607%E 060 


Oe 6, S0eet eee 
tees 0, 


Oe 6,07 Seecue 
tane 0. 


Ps o,GaS0tedn 


Ve 0,49eet eee 


Ve 6, sud0t ove 


Pa Hyd) bekoun 


ve 1.93008 065 
ue 8, 


US 1.9302 evs 
@ 


Us 1.SSenteus 
we &, 


US ,d50Sf 004 
ve @. 


ve 1,S900F oes 
os 0. 


UE 1ehet teeny 
HO be 


a 
q 
i 
’ 

# 
4 
4 
i 
‘ 


ccermt ome she Re oo we RA scene Tarai Pe oes 


Me ceaghht 


sees 


Sieteamstets - 


A - 


we ms eee Sot Ets a 


Sree le 
seit 
sorcit 
Srecit 
svecit 
orecie 
S¥eCie 
svecit 
Seecit 
aeecie 


Soecie 
seecit 
Ssvecit 
erit 
Srecit 
seecit 
seecie 
seecit 
seecte 
sPecit 


ove Cit 
orecle 
ecie 
svelte 
Swecte 
SPU St 
SeFcle 
srecie 
reecit 
omecle 


orecle 
. SUrULE 
seecte 
rece 
areal 
owrcite 
Shecte 
perce 
prec ae 
Serele 


reat 
wweule 
oveule 
yerdir 
bebtte 
Seow 
oreqbe 
wreble 
peedae 
Sebule 


horde 
set ye 
woetle 
awecte 
Settle 
pee te 
vetcle 
Serle 
der cde 
Seercult 


@ a8 1, eaegtess 
Has .1eO8e OOF 
wts etekueese 


tu 
coe 
Ld 
bad 

LT ad 
eu 
ne 
iy 
ve 


SAB PF eastvesny 
‘aS = 9 1e¥or cyt 
FIs .bOdH e0e 


fu a8 eee tert) 
has gidueteul 
PIs ,leusteve 


Lee 
Cre 
La 
ne 
ur 
mele 
1? 
’ 

"Ne 
a 


Bb ks teetdseous 
ad ,bseoeeus 
ws gtmSat ene 


Cy 
tne 
” 
te 


un 


2 


te az Peeeuge ony 
“AR ototabous 
Pie .tAtoe one 


'¢ 
a 


is AS tylawyereny 
ee gtdene ens 
PIs ghrude tae 


c: 
Lew 
. 
me 
aed 
wen 
Td 
‘ 
ue 


eS 8, 2270b008 


Puls L.cezdseeds te 2,7e90Kee) 


CELTS .0050R468 He =. ORE0ESE? 
hmte Jeleteos Sus 6, 
Cs 5 3000e 008 as e2iadiete 
Ce .e0e0toes is 42700806) 
Ce .03¢0toes BB o1SOCECSS 
Ca .deutteus ae 4 2dOE eeu 
Cs dover oed aS ,9053EoRe 
Ca {30008 006 as .2evseoue 
Ce Ke =,207SECCO 
Ce BS =. Sendtoes 
Cs ,drereove ta 4 3alitote 
Ce SNH OCS ts ,Sasetons 
HS OeBPTAE ELH Pris §,aagetoes 
ORES eeunnfest me ,aonk end 
Hie f71SPFOOS BKB Uy 


Cae gd0nte eve aa di2Stoug 
CE .e000r COS as .2?eetoes 
Ce .0300Re03 BS 4 t20MEods 
Ce .eacutoos as , fbb eee 
CB Phong ong as 4 $033to02 
ts ,t9hrecon 3B .eUGSE COD 
Ce .Meurt ety 3s eOTSE 008 
Ce . HONE COS 3B, VeOREote 
ts -Seddt eos ne, 3a71te¢a 
C8 gthyut eds as 4 5eSabeee 
HS S,GghUt oOU Puls 3,e232tou1 
CELE HONE CHO mB eeueteur 
mHus .TUb2FCO3 SKB 0, 
Ce SUF OOH as .elertous 
Ce .huntbeus “xs yeluMreds 
Cs ehimeeous as 4 tee¥eous 
Ce aeanntent as 2hnokeon 
Ce .de0ur ove xs | s033€oC2 
C8 .1voeeun xs ,2u9SEeoo 
Ce hPutteny Be LP ITSE O00 
(ore 3s WnoBFeod 
Cs edrulteud Re og fusieoita 
LS athoneets 1S 426246 00® 
Ws w.0e2rt ety Onis 1, susstous 
UELTS ohOSSE ONO HE QRNOIE OUP 
CU Seed Oe ee 
Cs .S00ve¢0u as ,P1edtovu 
Ca .e0untev} ae 4 eTRGEeO1 
Ce ahSOuk eds ae 4 teedteoul 
Cs eairbeul ae ,23nvt eo 
CB gervib one 1s, buSSsE Cue 
Cs ,tOvut eve is .duvsFedu 
Cs .weret ewe ae 429TSE CON 
CB a sAONF ev xs, Ynebtooe 
Cs Seutbens as a taliteva 
C8 gh bGE bens as 4 SOSGE OND 
HE HeOhleteieu HHIs F,vosVeebs 
VELYS ehOSSEOOH RB Ohad Od 
Weis 29 Mens Seo, 
GEG Suen ont ae get eet eow 
C8 ohneteout ae jdluuteut 
Ce .bitutens aS 1eGHECOE 
Cs .eeucbeut Xe .esduteuy 
fs Peat ote ya, 4N3SSE one 
"S gh¥orbere a8 ,@uvsbeuy 
CE HPL be ery ae 4 fISt eu 
Ls afhutbeus 3B 4 Yneoteus 
Ls ,beutbeve Be Sts be ous 
CE gllobbeus As SdSat On 
ms Pedlvitere tiles 2, luurteny 
URLS OUNCE ONY HB ATVAL Oe 
wis shhh Sbovs Ses 0, 
fs oSuthory 4B shh Abb ern 
CB ,notuupeny sa ge Fab eoty 
Ca ahgertens az gs euSteoul 
Cs .cdurbets as ,esnvteny 
Ce ,entngene £8 s0SSb 002 
Ce .$90GE 400 aw ,2u93E400 
(se .atetbedo as ,29iSt sou 
Ce atartiecs XB 4Wnbhb ony 
(sz as 4 fasdtoua 
Ce ehbuorers AB 4 SA58F oU® 


Hs §=,38028005 


Pe 0.1930E C00 
Tame 6. 


Sumbuts .0ss1bets 


te 2. 73¢utees 
aie 3eseeees 


Pe BOSE CEO 
tans @, 


SumuUts .BeS3Eeee 


AUIS Ue 
aduls 6. 
adults 
evugis 
aul 

abouts 
avuls 
eovls 
aoots 6, 
abuts 0. 


Va 2,7320E90S 
his ,3uIse 


Ps ©,0030b+00 


SUMOUTS .1080ECOS 


aluts 0, 
avitls Oo, 
AUIS Ve 
avuie 
euuts 
wouts 0, 
OUTS Oy 
wuts O, 
aviis a, 
auuls O, 


Ve 2, 72p0tee3 
His .s61vEs0s 
Susuils 


MUUTE Ue 
wOUts U, 
ALTE Ue 
Puts 0, 
outs Of 
AUUTS Us 
avits Oo, 
PULTE Ve 
ENHTE Ve 
Hulse ty 


ts ¢,7190b003 
“le ,sbuurens 
Sumyuts 


ALUIS Ve 
AUUES Ve 
avuts 0, 
ALUIS O, 
Atte Oe 
Ate ve 
enuts O, 
eugts vo. 
avats 0, 
AuLTR Oy 


ts e,iievreud 
ots. $nndbeus 
Suis 


fulte 0, 
PO eee 
ALB Uy 
avpuls vu, 
wouls oO, 
avuls 0, 
routs oO, 
enna Uy 
aude Oy 
avuls ou, 


Reproduced from 
best ava 


Pa 7,8950b000 
thas 0, 


eheTureny 


vs 7,778uteno 
Tans 0, 


oS SG0E 004 


Va J ee sue oun 
tans ¥, 


othede ens 


tlable copy. 


Us 1.50906 ¢05 
oe 6, 


us 1.59058 405 
us 6, 


us 1,5401t 005 
us 0, 


Us 1,56208 005 
os uv, 


Us 1, S70Peous 
Qe 6, 


118 $5 /M98 OU4 
ust, 


Samide aE 


eR 


Abus OS 


a ne oe 


%» 


thd 


SO 1,8760E004 8 =—« Bas 7, 000CE COD «= Pato 2, 3eseRoe? 
Me 4133SECGS «=F LVS §=,6055E C80 ee FORCE CCD 
Pts ,bonsece? ume .setegees Ste 6, 


weet w cs .decekess 4 odtatecee 
SPECIE Cue Ce .0000toss o8P00E 00s 
wETIC on Ce ,o8eetoeS§ = Bs. 120K 00) 
srecie m2 Ce .2auesees ne 423008000 
SPeCit ue C8 .200t oe? ae 30SICO? 
eecte = mdy Ca .sveekeee = keg 

MRE ne Ce ,6200E 008 as 008 
precise a Ce ,7euegoee ae 


090008 
Ca ,Se0etese Be, IIb obs 
wettis a C8 .theetous ae .S4detote 


AS LohOUSeCes =—«-_- HBB, HEKROVS «=—s PS 2, $2000) 
the ,tdatcees §«CELVS 0050000 he =. ISRO CER 


vie .byeateta omboltoos Bae &, 

erecié Gu Ca. sever eae ae .dldeve 
Srecle Cul Ce .oueetees Be 27 eekwus 
ercie Ca .oSectoes He 4 1200S) 
SeCie na Ce .haeotoes ne 425008000 
MPECIE Un Ce .200es oe? ne 8030? 
Srese rat CE .teneteen§ = ka, de OF C80 
Seecte ae Ce .Beurteve Be. aUISE CON 
wecie vu 68 210008 eee Be. Ceeetets 
SrEcie U2 Ce .Seeetoee Se 4 3aTItooe 
eeecle oa Co ,§eetous Be Sedetote 


Us UV, DeatEony hea wm, 97208 000 Puge 2, 7udstens 
Men ,NsleCt § «UELTS HOS3ECOO He = ,OChe CO? 
ws gtdsveced mele cazet~ds S5s 6, 


ovecit: Cu Ce ,Seeerece ae .2t2deeee 
Seecis = bne C8 .o00oFoul AS = 270eEoO) 
Soecle om Ce 4h Shot ons Ws 6120 ons 
Seecig ma ce ,dacetess BB odbbCL CCR 
SPECIE =u Ce .2ee0tee? ue. 503Bce2 
SPECIE = ndu Co .seneteee Be, k090E 008 
awecie wa C8, R2ANE O08 Re =, 207SE00e 
svrcie ou Ls .?evurece mB, fenetees 
SeCie ue Ca .deurecse xe, Se71beee 
srecit on + C8 odbuekeOd Te ty 


Ba 1, 1200F GO) ne ©, S54eseoy «= PLD 2. 00SOEP OS 
mas .th7%ee) §=6—DELVS neSSeeee He 4 Os 208008 
bis ,Iateteud whus .o2778e0} Ste @, 


SPCCIF Gu C6 ,duutt ooo Ne .aladteous 
Becer Crd Ce .euunt os Be el uekoes 
secie f& Ce ent oes Teri 22) 
avrCia ne Co .daenrous ae 4ddeet eve 
Selie oe 68 adberbeou?g se, Wa sous 
arPecie meu Ce. h9ORR COO Lo) 
BPeCiE x? fe -GDO0E 006 Ne, dO7GE Cee 
seecit ov Ca ,1enteove ae =, ees ote 
SrHette te be -S00uE ous se, 847 1Cebe 
SPECIE A Ce NOME ody Se gvetetede 


48 Iotiwetees Wa 3,03 5et00) PUle 3, ceo"t~ds 
WAS ,feuded) OCRLYO OGO1ECCe NO OOP IE Cee 
Pls ,1a9shoad Wns 01336063 Ske 6, 

aeecle tu Ce A WHA Oy ke gdteseeuy 
Hee «Lue Ce 00000 00) ot Me a te | 
secit om Ce .0300f oes 5B 128% OG) 
BvFcit me Ca .gecntee) Be ,23806000 
Brie ute C8 denn eue se, tusstoue 
Seecte nev Gs. 19GUt e980 Ya ,db¥3e 008 
Soocle fe Ce ,egbte run se 4 e¥ISECUU 
areule ov Co ,PhOnF age 5a, MeoBEoge 
Seecie vue Co ,devokode we 8e7 1G 08 
SPecie «6 Co 2h1eet oes Be =, Sesegees 


BS buMset oul *O t,07/00beut Ome 3,d00Steu! 
Mas gleewer} = VELTO ,egbutece He ,B7OTE CE? 
via ta7uceud wee 070800} Ste 6, 


BPeCie yu Ce .Judur soo as 42taSeeue 
eerclé Ceig Ce ,b000F 06) Ba dl eegous 
arecit on CE 45005 005 BS ,82e%Gou) 
SrceCie nd Co ,dtgekoe) wa 42380008 
Srecit un C8 .ceev€oud as ,30d3ECUs 
BreCiR meu ce .1teee cee ke ,2003Geud 
Seecie ote Ce ,s20nt ete Re 4a¥TbECb0 
seeuit ou Ca ,Pnonteye Je =, Veest oud 
Seruit ud Ce ,dervreve te ,be?itous 
specie oa Co ,Mleogees ae beSetote 


te 3, 7eeegees PO 9,00808 000 
OTs .361egee3s tans 0, - 
Surwute .atdeteed 


ts 2,0080gees «= P82, 1 Sue’ 008 
te .80128es thas 0, 
aumbuto § ,230atees 


te 2.0710600d = Pu & 09S OE 108 
We = {87eeECeS §=— Sane 0, 
BUTS = bINaE COS 


ature @, 
avuls 6, 


fe 2,00008e0§ O. PiSet ove 
nts 3703k003 =6tdne 6, 
sumputa =, 35318 +08 


ts 2,0eeveces Oe e,ekeotove 
mle .$793h003 = teas 8, 
sumpuls =, degen ees 


ve ¢.$0setees 
we 6, 


VS 1,0001teos 
e 


18 Toot sliEced 
ws 0, 


VO § 62008005 
Gs 


Ve 1,0078800S 
ee 


Ya t,e705F e045 
we O. 


3 


a 


4 
” 
a 


+ FE aes 


A yetn ct 


his 38 


sets je ETE hele A 


ey 


| 
4 
| 
\ 
4 
q 
‘ 


2v 8B BUT Ieb ees We so 12OtECet Onis B.e207sees fe 2.8000ECOS PS 0, BNTET CO ye t.oFettees 
ats , Meet oe. 


fi 
£ “As .1eRECOd §=—UELIS 05H OOO He PROTEC? 003 6 taee @, 
* His stavites2 weda 8802 oe$ Sis 6, Sievuts .0599et00S 
+ 
* 
7 
a ‘ seecik ww Co .senegese se .dbavecee evol’ 6, 
a SPECIE Cue Ce 00008 <0) ue .2706kos epote @. 
Z srcie & Ca .0306toeS «= ke, 128CF OO) ITS 0, 
3 aeecit = né Ce .2asetert as 42580t eee euuts @, 
Serie om ts .2eeetee? as 2 S05Sou? NTS &, 
seecie nev Ce .89est cos ae ab0secee evuta 6, 
SPECIE 2 Ca a200teee Be p207Steed woete 6, 
SreCae u ba .Peeetoes a8 euute 6, 
SPECTR ve Cs .Seeeteoee ae 4 SaTIESe abuts 6. 
sercie = # Ce .8euFoes as =, SeSeteobe ads vd, 
at a8 Svduseeul ws 1,8050¢ 003 ents B,007eeess te 2,.d0s0nee3 Pe 9, Cpveseee US 1, 7e%et eed 
Pas .108uteOs «DELVE 00628000 me =. 12O1E0OS te ,3080bees shaw 0. ws 6 
1s tezeeee?2 «meue Seizteds S¥s 6, Suwmuts .67608 0S 
SeeCit bu Cs .30unt eee ae ALedheve ovules e. 
Srecke COP Cs" ,ceerteous Ma Se7eG¢eus = ‘anuTa Oe 
attic & Lt .es@vteous as 82% ous adits 0, 
weecit me 1: ,@eertees us ,2Se0Eeee AUUTS ve 
SPeuae be CB kovoeose a8 ,SuSsteodd ADUIS Uo 
seeCte ody Le vet eeu ae. 200Stoue eouts @. 
Shecte na Ce ,ddveteee ae odOISECe QUUTS 8, 
ertcis ¥ Ce .Pevetous te, 000etoes auls @, 
SrCle v2 Ce .Seeetoes us .3078Eo08 woots 6, 
Srecse a te .tseetees as SeS0Eode wpOTe 0. 
a ee AB p,veneeeds ms togelsteot Onis $,76/eko0s ta 2,55006003 Ps S,ebuet eee Ue 1, 2eletees 
" eas oSPHeOUt «= PELYS «oHOSRECUY MB 672838000 wis 37908003 «sane 0, ws 6, 
a Mis .1S0seoe2 «Anus edSudteus Sune 6, Sumputs , §10Tteed 
¥ : 
a sercie bu Ce .3nonteto = Ww Cas aSECGG = UTS 6, 
BerCie tee eoOuvtoop” xe ,ePOuEoot evuls @. 
ES svecte on eoS60ESG$ «= KB tZOVECKE «= DUT SW, 
wi. eecst me edaevecot us 4dd¥0Eou0 aodte 6, 
. sreecie un o2e0ut 02 Xe, 8033EC02 
a ‘ SPEC It meu eh@vdeeee ne =. 20938000 
SPECIE ne oS2U0F oH0 as 42075000 aouts 6, 
a specie ou e2b0Khooe = KB VRORECUS = ALUTS Ve 
% Shevie te Pe ed £8 4 3G71EC08 evUTS Ve 
4 Srecle A ol lv0t oes ne 4Sebebooe aDUIS 0, 
BS 23 as Fvoeut Cos he t.d07dteo! Puls B.0O7TECOS Fe 2.530800} Pah 800CKCOR Ye 1,7081E CRS 
x Mas GAHSOREOL «= UEL VS §=,nONNE ONO he ©, 7725EOR) te 3700003 «= Fane 0, Qs 0, 
4 Vis ghseurec? wris .Speutens Sxa 0, Sumviits § .$neseees 
s SPECIE CU Cs .3990E +00 us .2laStevd AUUTS 0. 
3 SHEC IE bhe LB .e400F oO) ae ,2704o04 ad0ts 0. 
5 sreuth of Ce ,otultens Ke, teaOFoot WUUTS 0. 
¢ seecde ne Cs pduuteets iT merry sal) KDITB 0. 
4 SPECTR bn Cs atovut en? ua, 40SSoue anole 
: Meeudt Neu Ce ,teuckeug aw .2093€000 abule b. 
‘a SPECIE NO Ce .e200e ooo ae =, &¥7SEs00 avuls 
¥ SPrgdek ov Cs .7800F ous is = Meeseoee wdUls 


x SPrLTE Le CB pSHloFove as sa7iCeue RUITTS Vo 
i Seec te oa Le .dtereeds Bo. SUSREOUO AUVUTE Ue 


do kB HeHevdLeUD he b.353eb401 Pols a,teasteds Is dg suduteus Ve v.Geduteed ve 1.b0772t00S 
4 sas bOUIECUL «—DELYS onOESEOMY We e,aSLIECN? = nTB  s0aRR OO) = Naas e 


Pls ,1939E0O2 Wels ,eBHabet} gus 0, ‘SUMLUTS §.000R00% 
5 
4 arcie tu ts gduerb oud ae oelddeouu anula u, 
: seeCle tue LB .munoreus ke ,27uatoul 
Y Shecit or Ce ,03008e0% AG ,1ke%Eous 
Shecte me Ce .esuveov) NS 425806000 
. wecte oun Ce .aaveteoce Ms, W380? 
; Prod LB .h@undeue BB .2u9SF eyo 
4 SPeu dee Ce aeOEF eno aa gdvsbb oun 
y arecit ov (se 7 Bbut ous £8 ,Vonktooe ; 
SyrCie ie Ce ,Snitteds as ,3elitege a 
Swecie a Ca .180bkeous sa, $0SeE=0e 1 
rh) ae VSI ESe oor mB Lee bat ent onis ©, s2nereut ts ¢,dututeed es o,0ugureod 1B JP Vouse ony q 
vas atovateut PELYS onset edu we ©, /A¥SE 002 wis .sederees faae vu, us 0. ; 
Mis clevuteud amis ,8ngutou} S18 Uy Sumpule ,edesteud ; 
: 
secee GY Ce Wout eve pe @taseeou = wttite oy ; 
seecht vue (8 hump ent as ,2fnepeont 
5 seecit oe te .eseubens ae ,1eeveou) 4 
owtule one La ,dauoreut Bad tnveouy 
a SPEC LE um CE odeurden? as ,sussteud ' 
orec ie re Ca bob ere ae ,d0Dst eon evils 0, 
q Ole ne Ca cdetuteus «= as 20ISteuR «= aut S ue 
i aretat LB oFHOMECNa = BB VORNECUS =— aUUTE Uy ; 
; sreult ote Le ,Sevnteny ss ,Su?izede eVUIE Uy ; 
% sveclg & Cs .d100to0% Ke .5856boue adits @, 


B-5 


oe Deemed dant eh tenn Sangre inn ad Ca tae LEIS Sa EES ATES 
errant materi iS ene 


one 8 ,tkevi. 
Ste B® oibevi, 
@ 000s 

© .350E- 00, 

8 Se one, 

@ o25ebeee, 

S odtowls 

S .Sederuds, 

& .ddsbees, 

SB $0.0. 

8 ede 

8 & 

St 

& 16> ‘ 
MOOS, c2EOSe oSEPODs cOEOS, Oc0e GeOe Oe Code 
SB OeWe Oe Ves Gee De0e beOe 0.0 O60, 


+ wt BGeWe Vee Goto O40> Cede Ces Vote Oe 
198 BS OeOe Vebe Vebe Vobe Ober OeOr Coby ede 
Cd | BS CeBe HO Dede Oey Coe Oey Oo0s 0,0, 


“Firgomine ww 


e 4 

PONS IAEAM VRLUCITY 1.G70F 00S 1,300ECOS «i, B80keeS «= 1, 2eOE COS 
CROSRCSTIUEAM VELUCITY —d,7008002 <S.S00E CUS 7, 0808008 91,0108 003 j 
PARTICLE TemPERATUNE 2.H1SE00$ F,B03E0OS 2.023003 2, 97REC0S 

PARTICLe VEensiiy eeTVRoUS «8. 6F0teEe =, UaORCUs «=, SWE oud 

STREAM INE YU, a 
Dr Si véLucity EeMTORCUD S.F0F ond teSdOkoud ft ,ddetoud 

CMUSSeST NLA” VELULITY 92, 90b002 ed, C2EF 002 98,14608002 1.070603 

‘PARTICLE Tee Etat UME 2GIHOes FwHKECeS 2.055003  2,073E003 


[PARTICLE UENS. IT 2ePVOFeNS 1,030F ODA = BREF HOE «= SEE MOD 
catwe Am int Ue s 

*pOmnstnc am veuity 1,GBUECS §=1,S¢OROUS «= SAOECUS «= Sa SOR OOD 
‘CWUSS~8 Tuy Are VELUCITY eS,UutOOd ebe2¥0RCO2 HO,000G002 91, 10000) 
\OAATICLE Tammewatuwe §=—«- 2,8138003 §—2,0958003 «2.015803 «2, 078003 
«pant ice Orasiiy 2,790E9OS «1 3ROEreR = BREEN =, FANE OOS 

STREAMLINE NUL « ; ‘ 

DUANE IMEAM VELICTIY He OBUFCUS 1,80CECOS Lo BaOECOS 1.dheCeus 


CHysSogiiean VELL ITY oF, SLOG M2 CO.7N0E COE o%.220beud 91, 210800d 
PAUTILLE lePvedatune PoAlStoes 2.093E 003 2983 ous 2.97 S005 


MARTICLE Uswotty SePOOLeUS 4 ,850FSOe = APOETOR «= EOE OOS 
STE AME Ihe ws, s 
OUANStNkar vELUC ETS s-MOteND §=1,00N800S = SAOECOS «= 1 HOE OOD 


CHNSdeRIHERP VELUCIIYD ef, d5uk002 O7,999E 002 OF, ROE CUZ 91,3008005 
PAHTILLE Thowgwa tine 2eO1SCOS «2, AGIECOD «2.OLEECGS «=. TE COS 
PONTICLE Vensiir CoWPehoeS seASeEwee § §,820ECRe § =, 9408S 
SINEAML INE WO. e 

OOMASIAEAM veLUciit OBORCOS 8 -000ECEH = SZOECHS «=—«-_- 8 THOT CNS 
CHOSSCSIWEAM VELIICITY oJ,APAbCO® O7,790E0F2 oh .UOOECUS 91,4008 005 
OAWIICLE Feaen atime THEM COS FOOSEOUd 2.09EouUs 2, 0738008 
MAWLECLE UdnsttY eeMWECHH 1, 430bOOR § 1,820¢-08 $5006 eed 


SIMUL AM Ime AU, ’ 

OOansingan veLUcsit eOCOECOS seAOECOS 8. S20ECUS «= BOT ONS 

CaOEseS ita VELUCIIY OB0OECO? 01, 1808003 Oh, 510K CO) 
FBASTOOE «2.913003 «8 FSR OS 


PARTICLE IemmEwa Time 2. 
Pomticee vewsite AMVC VS FeASCROUE 1,020EHOe 8, S00ke8D 


STH AML Ine NU, a 
1OOUECUd be BIORCOS ne STOECOS «= BOON OOS 


OUansing de veruciir 
CaUsdeBintam veLUCIIY OO 1GOECER O1.2008003 of ,0906003 


eee cea) ser nan bate ae 


PAWtsCLe MevmeRatunt a 3B -eBORECND =. 20 913ECOS §=—- 2 HOE COS 
OONTICLE UENELTY 2, 79VEDS =H AAFOECOR «=, B2ET COG =, SUE OOS 
SIME AML ENE Wile , 

QVedating am er Lciiy PeSUKEONS = LGIGECUS = ELTSOR COD = RODE COS 


CHUGSH HWE AM VELOLITY SU,G7UECO? O2,080ECOS o1eS00E0OS of, 0208005 
PARTICLE teewekatune 2eHUSECOS «2.003003 8,013E003 2. 90REC0S 


et ea Ook 


PARTICLE LENSITY 2.7O0RUS §«1e4N0EMOe = SROECUS «= FAECES 
OTHE APL ING ai, Wv : 
Ouadelbe ds veuUci ry 15008005 1,620600S beS80RO9S fedoueess f 


CHLIBRCS IME DH VELILTIY eS,S2Ub O02 ef eS2ueeOS ef ,d00E0US od, ugue ees 
PARTICLE lemegnatune 2cF¥St00d «8 07ER00S «= VORECHS §«=—-20 9538-003 
PARTICLE OENSTIY C,TUOECOS =e ATEROR «= OTE TOR «=H FONE MOS 


4 
4 
1 
i 


ce dite SAS tile 


a 


gia ANY 


Sy 


eee 


' DWE INE ay, st 
OOebSIWEA vELOLIIY 
CausSestecam vELUCES: 
Panticce temPendtuce 
OaatiCce vensity 
SiMe dma Ink td. iz 
OGangtacam veLucity 
CHLSS@Stat an VELUCTIT 
PAMTICLE TenwEndtune 
PantiCLe veowstty 
Stee ams sae nti, 43 
OOens tHe an veLUCity 
Coussestntaw veucity 
PARTICLE Vernet watioe 
Pautitce Uenaite 
Stet emetag tar, ny 
OUehStatam veLuLity 
CHOSS>SINEA™ VELUC ETT 
Pauticre lemmewatuné 
PARTICLE OtnSity 
SIME AML INE RU, ts 
Date am verses 
SOUSSOS INCAS vELUCETY 
JOARTECLE TemPEMas UME 
PamTicce Density 
[STREAM INE WO, te 
WUenStetde vELUCETY 
;CROSSOSIMEAY VELULITY 
PAMEICER 1EMPERS I Une 
MARTICLE VEnS)tt 
{SUMEAMLINE wt, Ww 
OUahStagar vVELUCETY 
CAUSSeSI HE aM VeLULEt? 
PANTELEE Ferre fit 
Pawttere obustie 
BTMEAML INE mel, le 
OUAhSIRE AY vVeLUC EE? 
CeuBresingas veuc ate 
PARTILLE Teme na tint 
MARTCRe venwsity 
Stweamine md, Ww 
ODANETHEAM vELUCITY 
CHOSSOSTHE Am VELILITY 
PARTICLE Teme tal UKE 
PAnbACLe Vensity 

Stee ari tnt nt, ev 
DUenSIHEA™ veLUCi IY 
CROSSCSINEAM VELUICITY 
PAMTECLE TemPenatune 
PAwTSChe vevelty 

SIWR AML be Meer, ca) 
Duarstmeam ve LCi ty 
CHOSDeSTHP AM VELULITY 
PAWTECLE Ieevewatune 
PARTICLE Otwalt 

Stmk O4L Lk td, ee 
Dusratneah veLuchie 
CuUrSesingar verity 
PARTICLE Te WPENaT URE 
Par iCek VENSITY 
Stheamrint Nu, r?) 
OUanStweam we rol ety 
Caldsestng am vera sll 
OwPElte Trsertatiane 
ParfiChe vengttr 
Statamctve ww, ee 
VUANSE MRAM src TY 
CRUSdSeSImEAY Sette 
Marbgece bee oearttunt 
PamTECet ve ovate 

Vee AsLINe wu, e 
DUANSIRE ar veLUCI TY 
CuUsvedsiongan weeitaty 
Paws cde te bi watuwe 
PARDELLE ver SUEY 

uw oo 


1 

eu " 
Son " 
a Aa] uP 
S tv ‘ 
Ce) rr) 
Torate me 
eo we 
aon ue 
tet ee 


.310E 905 
oo. 21 ok coz 
20 203E08S 
2eteetees 


tS2edees 
of 000g ea? 
2.703 005 
2, F00e 065 


2. FO9F +03 
foI Mt ed 


tedvertees 
0,030 002 
2. Toes 
3. 190¢ 08S 


1, SSet evs 
mb epavteus 
Ze PoS ees 
2. Feb 00S 


§Seetess 
ot, onb€ 003 
2.7aous 
ee PSU ods 


1, S806 005 
o4¢,weuteus 
eo faseeus 
ee FMetouds 


1, S906 05 
24,91 OOS 
eoTpstous 
27408 00S 


ite 
LJ 
we 
hau 
Cig 
me 
Lrg 
a) 
is 
ale 


$.030€ 005 
1.2006 063 
2.073003 
$.0508eee 


8.6308 005 
3 ,0a0t 003 
2,873E003 
1.63etcee 


1,0008ee8 
©1 0008903 
2,073 003 
194506 ove 


b.edeneed 
9200208603 
2.699€003 
1,0$0¢ 0e8 


t,anee ood 
oe Seeteus 
20 BSH 003 
1-0Setooe 


1,870b00% 
93, s00ree3 
e0esteus 
§.aS0tooe 


$.49at ood 
oy.vase ros 
2engteeos 
t-aScteoe 


$9006 005 
#100206 000 


2o7O3t 003 , 


1.440¢ one 


awcrrrteae res 


003008005 
ote P ECCS 
2.003 065 
tOdtkoee 


$o38egees 
93. 086E 003 
2,003EC03 
Le ddegoes 


$0300 ees 
02. d1dEees 
2, G09E 003 
teBevtove 


tedovteous 
oR, 1I0ECOS 
2089Gee3 
t.82eg ete 


fodlotees 
o3,400tee3 
200898003 
tedaebeoee 


$e boegees 
a4. ete ees 
2.0730 003 
t.B20G ees 


1.S¢e¢ 605 
bo 920E us 
2.075008 
Leddutedea 


$,400t ees 
o}, Sect ead 
2o0SSE005 
1.020keue 


&.270€ ees 
o2.2Tet ees 
2.9398 005 
S.d008eeS 


8e8P0ECOS 
22,0006 ¢03 
20483E 003 
d.d80Eoes 


1 -B00ECeS 
o3,030E 003 
2.0098 003 
S.deveees 


$2008 005 
$0306 003 
2e933E00S 
9.5408 00% 


todvecees 
-8.Seerees 
~ 20958003 
$,5408oed 


$300 00S 
oo, 00veeos 
209250003 
S$. 900EoRS 


1.StCEoes 
bese? bpd $ 

ae 9230005 
S.daNE COS 


1, S20€eus 
od ,e20beoe 
2,905t¢05 
S.S00ke0S 


oh OFoqueu,y 
Pe td hd eo 
odVROEG0U.0 
o S0boOHCULY 
eS¥Eri2ou,l 


Maia tee Se saad ae abe ait oP nel tor barber std 


even 


ve 


082600 002 


as 2,2boeeest 


ass 
vis 


ob ssmoet 
eveuteeet 


Cu 
te 


« 
“e 
ue 
wev 
ne 
y 
le 
a 


BB 2,200eteud 


was 
vts 


ot FT ieees 
oVdlorees 


ae @,27¥oteut 


was 
ets 


srecit 
arecie 
seecik 
evecie 
avecic 
SPECIt 
aeccle 
seecie 
SvPUit 
Sorcie 


ot7Seeet 
oPZoet oe! 


tu 
tue 


Oy] 
UT 
un 
nal 
na 
uv 
ve 
a 


aa ded? Sues 


“aa 
Pts 


Srécte 
Sect 
avecie 
avecig 
SPECI’ 
SVECIE 
avecie 
svecit 
SPECIE 
Serle 


oh PSE O85 
o%e8 $eeos 


AB d,ceule 0} 


as 
vis 


Srecit 
SHLCUe 
srectt 
aeecit 
swecie 
srocle 
Avec le 
SrEC It 
De te 
sercie 


ott TSéeo 
ePe tubes 


as dgeduncery 


was 
vis 


aerute 
dvecee 
seecit 
serch: 
Seetae 
seecit 
Sercie 
oreeit 
Srecit 
seecit 


eh f sates 
Vb eeu 


ee 6, Outae @. 
As HH deearehs Puts 2.240itot? 
URLS ,S22dECVO HS ©, 7110ECO? 
tye 48077063 Ske {8621ECC? 
Ca .2veteeue as .2t0egece 
C2 .nesreoet se .d¥OSEoO 
Ca .3123Fo03 as ,cedtteed 
Ce .26ent-e1 te =2050e 
Ca .0d0et~es As _S90SEoOS 
Ce .980SEe0e Ne 4 2090E 008 
Ce ,ezeete Ns .29ORE 088 
Cs .827Tto as = 1S9UECOe 
Ce .beatone ae FI23Eoud 
te ,sivegoss as 4S477E0 
Re 1,003%000 Puts @,cehetoe? 
OELTS R20yE 000 we ©, 73528002 
MUS .SNMSHOOS Sas .feduered 
Ce .etose cee Be 4 2teneeet 
Ca 20ST oO) xe 4 8070E 00) 
Ce. Bh8SEe0s Rs ge! 202 
Ce ,2hertoos us .d437b eeu 
Ce Bact onys ne ,vesseces 
Ce .tever cue Beg hhVOECS 
Ce 4200600 us .20et eee 
Ce .t28et ove Rs tS07ESCe 
Cs that oos 18 9 OB19F SOS 
Ce atte eos as 9870E~00 
ha 2,0d¥Sten9 PMs &.119Ub—be 
OELYS .PIORF COD He ©, 7272008 
AHOe .308eECOS Sue .10zekee2 
Cs .a¥nit 00 xe .d@hloGeve 
ts ,oddeten) ae odMedeeus 
Ce .30S8t 003 as ,ettete 
Ce .28a5bo01 ke 84388000 
Ce .a5e2e 003 Xe, 10038202 
Ce .1e9AF eu Re 9dlvjeedo 
Ce ag eco xs 2veeteun 
Ce .tdubeous as ,tagheooa 
Ce f1YME COs Re 4 elOvEcus 
Ce othturedy is .b42nEoue 
we 3,2e7%areue Pris O,ustagoud 
VELTS ehtRUL OUD mB o,7eu2bou? 
Whus .SO7TSECOS SXB .s01Seud 
C8 ed4ervFedy aS 42thieeno 
CB ehStTben} ae 4 @4Sutout 
Ce edunnteds Ke ,el0ated2 
C8 ,aaserodt aS 9 d4S0E 900 
Ce .0b0hto03 ta .t0lieobe 
Ce .suvareco we 428016000 
Ce ,e2ote enn ks 424008000 
La ,1saebous as gteSiFouve 
4S olbbctons Ae ge0Steous 
Ce 4800FeUS 4B 4 da7UECVe 
HS &,On97t O00 PHIs 3,12TeEoot 
OFLYS ,ATASE 200 HB ©, 72H0E Cue 
wHUS gJuerspenS Exe ,tulleeud 
Cs .2972¢¢00 BB 81148000 
Ca ,eahetens as 4 &918tous 
Ce .30387¢ 003 KB ,eUS0ECUe 
GB oduditeant as eudstt eve 
Ca nrs3sboos ae 4 Ueseoud 
Le ,14Gltooe as adldutouy 
Ca .4260te00 Re 4 aFHHE 000 
Ce .tserbeds ka gL 02E008 
Ce ,te7eFede as ,neheto0S 
Ce .1tb0hend re gD KEeve 
eo 6 ,8V0KE SOU vnie t,sedteout 
NELYS ,OI7TOECOU Ke ©, 73Sete08 
hmus .S0STEeuS xe ,1019b002 
Le ,eb7uteve ae geltorede 
Ce ,eeuve yt ae dWwitevt 
Ce ,3g8tCoS ae nOheeode 
Ce ,Paurtedt Ba =, 2a30E ove 
Cs .B74Steus te hueerene 
Cs ,s4oteeny re ,2lveteuw 
ta ,uguuteuy Is ,2vangeun 
Ce ,Lds7Atous ae atitegedu 
te ,bustteoe ue 4 059St00S 
Ca .tiveteus rs ,5e78Eo0e 


Ts 2,es0S8e05 
oOtiteees 


ots 


te 2,en20eeed 
o8107E003 


ts 
Suewuts 


Te 2,0$3etces 
Py S23) 


“te 
Sunduts 


ents 
avols 
wooTs 
outa 
wonts 
abuts 
wits 
evuts 
avute 
abuts 


ts 2,uyedteo3 
A sry 


wis 
SUNOUTS 


wouts 
avufs 
avuts 
auuis 
avots 
wouts 


abuts 0 


ays 
aulits 


wouls @ 


Te 2.45306 003 
oG405E003 


“Ts 
SuMvUTS 


n00le 
advise 
avute 
aulile 
ants 
avuls 
eoute 
eod0ls 
edits 
auits 


te 2.asa7eeus 
00955003 


“Is 
aumouts 


auvle 
auvuls 
evuts 
avols 
auits 
ayols 
eulls 
avute 
avis 
adits 


08 38egeo3 


081M eee 
0081 $08 
ass 


ED 


OOF ke et ed 
030008008 
ot237E 082 


ot hPOtoed 
ol 2lotoe? 
98180too2 


o11223e0ee 
oS0s2toee 
Ve 


oS PCE COE 


o.3useroes 
S060 006 
9002085 

~.0S12E003 

©, 9oSetoeS 
oOad0E~es 
e. 

© ,9357toeS 
e2OnOFooe 

a 


oft O1E0e8 


200008 
e2OVOE Us 
2991 7F 088 

w.02eatods 

o0577t003 
0% 72003 


2, 898SE 205 
e2Slotove 


oSeSht 008 


©. 16S0Eoee 
otaSetoes 
050384003 
SIMO CHS 
oobi Parous 
ode? Meus 


Ve 
#,47dt00S 
of 5 
U, 


o8eldEeoe 


o.bse7t oS 
areoedy 

0 3OO1EWOS 
o. bezseces 
oe SOHOFony 

eS9tebeos 
Ve 
esdiboteus 

600SbC0? 
Ve 


Wu 5088 
vane @, 


Oa S.813dE00@ 
tans 9. 


Pe 3, stendese 
taee 6, 


Got ote 
ao ge 


ve 3, 1087 obo 
Tans vu, 


Pe 5.09756 00 
thas 0, 


Va S.08eltene 
taas 6. 


us 2,00%68 00s 
ws 


bs 2,0izitees 
oe 0, 


Ue 2,0180C 005 
eo 8 


Ve 2.0ud6 069 
use 


UB 2.01008 005 
ve 6, 


US gevUFIE els 
Qs 6, 


Fee 


Sh os 


er 


. 
+% 


ear a ll 


Sra aa rn oe 


‘3 


Be 2.2027t008 
mas 200 sseeet 
Pls ,Viddeee! 


Seecie Co 
Welle = Cie 


é 
"Eo ESE? 


* ne 2,2duerees 


man J APT It 0O8 
Pies ,Cevteces 


co 
cua 
C) 
me 

SeCle = On 
au 
na 
u 
v2 
a 


ae 2.2030b00) 
Ae 48 7O0ECOS 
Pls .eulacees 


Seecie Cy 


KS 2,19Set 001 
“as 4b lbeteu) 
wis ,Beseteos 


specie Cu 


SP;CTE ue 
Srecit 8 


£8 2.17SSEeu) 
Mas gL PYOF ON, 
Vis ,bnsereny 


Secs by 
SeeCie 602 
avrcik om 
SPtLHe me 
s¥eCie oun 
SetLit: reu 
S*eCie ne 
SeeCie ov 
sreute ue 
seecte #8 


as @,tssdeeul 
whan ghiveceuy 
HIS MAD sbony 


seeCie cu 
SPECIP tue 
srecit on 

SPecit me 
SPECIE in 
SeeChe seu 
Serie ote 
shecie ou 

SheCit ve 
seecit «4 


fe STOCIECEO «= PES LSEZNECEL Ts 2,0517 003 
ELIS OL7SECES HO ©, 7985E C82 


wpe = 430akeoos 


CB .dOTOE CCS 
Cs O37 o8! 


CS htesteus 
te o.d353008 


COLLIS = .HIMIROC He ©, 7957E002 
te «= S02sko0s Ske 


Ce .2erercee 


Ce st3e7Eeoe 
Ca {teterote 
Ca .8seerees 


as 7, 38zeeoe 


VELYS .O1OSECRO Me ©, 77008002 
ands 


02000 003 


Ce .a97Ste00 
Ce 2058 61 
Ce .a0adEoes 
Ce pdeadtwot 
Ca AGERE OS 
Ce .1948E 000 
CE .W200F 000 
Ce 499308008 
Co .8OsaFoeS 
Ca .t100kous 


We B.12v7b eNO 


VELTS ,H221F 40 mB ©,6071E 402 


Wnye 62%e7boos 


Ca 4207S +00 
Ce .H3@ob out 
Ca ,2@oureus 
Ct Mey ad 
Ce 683706003 
Ce 4 1908E +00 
Ca ,aPuCk +00 
48 otehereda 
Ls .¥eunt~ds 
Ce dsudeeys 


WE 0,%2138t 000 


Ce ke7Se enn 
C8 .SOrF od) 
Cs .7hROten} 
Ca .2eaches 
Ca .heTar ous 
Ce ,uacevu 
2h2UOR CUR 
eleoitode 
Ca ,Gageeens 
Ce tb hkeus 


wey Pye0P ent 


Ce ,e@v?ubeuy 
260261 


(s cai aseeos 
eo 1 Qhat tou 


oldeatone 
Cs .O1Steeos 
Ca 4tivokeos 


Os 3.07288 000 
Oto ,8000EC8S «= fans 6, 
Bro .1013Eete «sunbute .0301fees 
se 281 7Eeee HOUTS ©, 1200 O05 
NS .2AOMw~os ebuts ,t20vEres 
18 ,d90SE~O2 BUTS 2150 o0s 
RO 4 2629E 000 avute 21926003 
me 10388 wOOTs ©2171 7E90s 
BS .2107ESOO wOOTe 2171E03- 
ae, 000 votes 6, 
as 3 ?a2tuee OUUIS ©,00302°R) 
AS ,e405F 06S auUTS |, 2197F <0? ” 
as .4aTeE~te ovis @, 
Outs 1.8107Eoe) fe 2,eanetees Pe 1.00006000 
tte =,e8ekees «= tans 6, 
efOtteeee «sumvUTs §.62S3E000 
we 12197808 adUTs ©,15106°0S 
Me 4 2087EoOs MOTs .8510E08S 
ae .S920Loe2 aduTe 21098003 
Ya 4 2adOteee PDUIS ©.2169%003 
Xs ,te1eode MOUTS @,220SECRS 
te .2uereced eects eens? 
us, TT) epote 
ae = tVe2eees nucte oc007TEbs 
as ,nS23b00S abuts o,194 76-07 
ae nalhEoue evots 6, 
Pulte 2,¢ez0foe) Te 2,0008tee3 Pa 3.0177E000 
Hle ,a0eSéeo3) = TAws 6, 
SKS ,8008ECO2 «BUMOUTS 16008905 
re 42t37eou AVIS ©, 2208E ous 
ue. 2003EeR) WUTS §.2208Ee0d 
ne .S808to82 B0UTs ,2634E~e3 
48 92830000 ADOTS ©, 242TE~ES 
aS ,8002be0? AVOTS ©, 206oto0t 
as gdlutteue aVits ,2e3ebeds 
as , 20888400 avUTS vu. 
RB .LeSTECOS RUUTS ©, SO3OE ote 
Be O163EC8S ADOTS ©, 130CE 07 
ka ,5078E000 RUUT® 0, 
Pulse 2.4153EOE Te o.4hO5t 00S We 2. 9801Esee 
Wis ,40dnEso3s «tans @, 
Sze ,1udsteod = «SumuTs .1270bees 
Ma ,atiteeon LUTE ©,2630E 205 
XS, kOVHECYY ours 6dus0F eas 
te obPieteve nuts .2eotteous 
as daditovy AUUTS 9.2uS7E0OS 
28, 9U0M COS LUTE @ede7 bods 
KS eU07e Oy aouls 463883 
Bs ,2@98RE000 wHuTs 0 
aS ISON OOe aunts aohtode 
xe SVALE OOS auits ye@rea? 
ae .Sa79Eede avuls vO, 


vHle 2,509%beu) 
MELYS ,AgNNP eNO we @  A2STECO2 
MMUS ee2Ieous BEB 


oi Ouieeo? 


edt i tkedo 
o8897G 00) 
eSTSAR CN? 
28S Fenn 
oVOCieous 
ofhulgedy 
o240ete00 
oS STOR CUG 
ebATYE COS 
eSOTRE CdS 


vala 8, H2niteous 
MELVS UePTreud mE eo Bavdte'e 
wHiS «= ,2HRIFeus SUS 


eWHo4t 204 


s2tlegona 
ePAVRQ oy 
ebO71€002 
ed432F ogo 
eFrdet ous 
setuleoun 
2VnGe oul 
st82ebone 
oS002bo0S 
2de?Weove 


ts 2,0se5teos 
nis ,60SAE 003 


SuMoEts 


wdUt> 
edule 
wiuts 
evuts 
nuuts 
avule 
abuts 
aunts 
abuts 
avis 


is ¢,edtofeus 
wis ,wuseteos 


Seyats 


ruts 
avuts 
agnts 
euute 
atuite 
nts 


Pe 2,0300k+ee 
Tans 6, 
ob S08 005 


~22601Ke0S 

a01E0S 
obSO2b08s 
wo2Sb Sous 
eodsletous 

od SS0E COS 
oo 
moO TI bede 
o, 962 St eue 
We 


Ma 2, HNSSkeuN 
Vane 0, 
oheddeeos 


eeteudgous 
oddMdtooy 
e@20ieoos 
@o223etoes 
o@@Slbo0§ 


odeSteons 


BY 


ye 2,00028¢05 
ee 


VS e.OdtIE Ces 
os 6, 


Ve g.0127800S 
Oe 


We 2.09078 00% 
oe 


ee 2, 0207E 005 
Os 6. 


UB ¢.VPNAE ONS 
use, 


i 
: 
i 
4 
1 


$e a8 eelenveers TE HLdESt ett wels S,leadout ts guest ess VS 2.825% 000 8 2,0319tou4 
“AS .1OUSOLE «—RLES Ate deer MS o,fbtet ee? Wis ,eeseteus «Thess o, ae a, 
Ole arsteeus ame gacherS Sas ,bUbeeO) «suments ,e8zegees 


Seem te Cae FORA ene Bs 42hteteoe arlits ©, 2$Setot§ 
weit ue Ce ,eauSt oe) aS .e¥COReuS wits § .25%0t05 
s’etle on Le edhSicgens ae 4 Seletousg wuts 

arecie ne CO .gesiteoes o2eddt tus evute 

sercle te Ce .Mhleteed 09500005 euuts 

Ue Bey t8 ~eratece eft ertcoy abuls 

SPeChe se Cz e200 cur eo LYNG C08 ants 

0 Le .derutore ebsdateue ATS eo hadgitevs 
wecae ne Ce .¥9Setous Pel bo t PUES ©, /oUSE oes 
artule oe £8 bs Cur ens 88 .Su7Mtove avuls &, 


{> = as eotvdverus Me L.UStiters Pris $.68376008 Is ¢.azancees Pa 2, Feat eey Te 
whe LIME ee eus VOLPE eHOOTEOUO ME ©, F9UYEONe wis ,eCnutees thas @, use, 
whe NSvdr end — kmE sePdibouds SBS oamtbeod Staut]e eesoteady 


Seade tn Cs .d4ieteuc we hi teteee e0uls <.22638e0$ 
SeFele tue Cs ,oattteet £8 ,200ot0us eduts .250oey 
Sercit on CB .2nzttons aa .Se21Eou? ants .22edtees 
Sel it me LB pddegten) as gdattteuy 

Sercie tm $B POT ot3 AB 4Veeateous 

Seedse mde C8 1% Stour 48 adlvltevy 

arelse ne Ce .Meuvteuy AS PCO CCC 

Semele Ce .fivseowcg BS. baMStoua 

wel ve C8 4 ehsurery As SoM 064 ALIS ©, e07uteyn 
Seefie oa Ca thenpeng 4B {dal dtoue aveits ¥, 


so aw eoulede owl mB Leelee vient ris Seimveteut 8 ¢.ehterous PS Peoseuteuy UB P.O%eerouy 
sas derieeus «VEL TS HSak oul HS ©, OyJetee2 ates Teoeds tans o, ve a, 
Te .ANONEEYE «HR QehOrtous SKS .VI9SECOD «SUMLils .2702b00% 


Seocle out (8 .POTe oun ae dt iSbove OUTS ©, 2oeut ons 
seeede lee CB peusotons As AE Stout evita ,dsauteus 
: Seecle bh CE 27 uStends Be 4 S402 BDUIS =. 223SECeS 
areChe He Ce .eduttert as 4 daSetedy AVUTS o.2dddtod} 
SPC ge tum 1S gPasPbonys BS ANZ aR ONY CHES eeeesEony 
weecie met C8 NGusrene a8 edlultevy etuts ,dedoteus 
y eLEe ve LB yweMiteny as 2teveete eulia u, 
precise ou 48 ebititeds BB g1I1Vdg oun elUls ©,65¥0t coe 
Swocie te CB Mavereus aa SeSSEous SLUTS ©,0030E°08 
srecit o# Ca yi targens AB SUhUF OLS aOls 0, 


$7 a8 d,ubreony mB Leeviereut brie a, tavurout Ts ¢,s¥safegd Va dodtlOteus UB P,0maleony 
YAR gthusteut VEL US hirer ert mE etrastends wis ,ausetsus Ves uv, us 0, 
VES ashe treed mms ehurhers Sas VITOR Or) «= SuMtEtS 6 StSeE COS 


Sercae te CE eit eun as ,20iShevy ele egestuteus 
route Cae CB puget eut BB gevdsreus tile 42$tereds 
Dee uin om C8 oertet ens as ,S20eteur rHyle ebaltees 
aercit me Ce .euadt ous Ne 2a tet euu .2usetoes 
aetcie un Ce .Miedeoos 8 ,MSSUF COS *o205S6 003 
sre Che mee Ce .bMGSE eu" a8 42lultouy 02039003 
Shee gh owe i Meee a BS A YANE ON0 CRUE Sry 

swrele 48 wh eb ony Ke ytessfeoue fue ©, fesbreun 
seecit te C8 aTaleteusS AB yeryoreuds abies eo teltoun 
peecit 8 CE ghhoteens aR SUBIR ete amulet, 

ta a8 eyutivarenl me 1d7eseert nia a.9/9G C0) Te 2, 54e3E O03 Ve 2.5071 ou0 UB ¢,b1bateus 

Aa gtaeebeut PELVIS Witerent me esttEeens mis ,wudseteut thas vu, uso, 
HIS Gib srend Wee ge Gtubenrs Sas ,Onyvreout Stocels ,ssorronsy 
SPecde cee Ca .e8Tibeuy me 42@histere ALUTS ©,1S 7B 00s 
serCit lite Ce ,eeSsteus Se .PVeetoul ADOIS ,ASTOECNS 
Ue Ce .edSerdens 48 NV ebene ws falebegs 
nee Cie ove L3 ,eaurdert it Mey Ba VUNIS a, felegeny 
wetye ove CE gm Parporys AB gt¥ytpens FUDTS eohutategs 
SPecit men Ce ,§8utheoun Be eld? euy obits .bategegs 
seecie Ne C&S .derbbeudg WS 4eVVUE Cul avuts o, 

seeudt ou CB WOM ECDs AS ,1isobete ALUIB ©,4900E C08 


nerd ge oe Ce amPipengy BS gdeettors Mille e,@L7iFeon 
seed ye & C8 gbilvrens 58 {Sat lbeon emda, 


1% aw 1, YENWr ony we t.adSd10 00) bala a,fesitods fee. siidseens ¥a 2.duteteue 8 Set stareus 
wae hviveous UELIE ,ONhOE E00 me eo, javetens wis ,s$¥autens fans &, WB ¥, 
Mla ,fucehovd = amus .eSeiteuS Se busted) SuMduls ,Segetens 


Seecue oe Ce chest ouy be glONE eoy 1S ©, Sse¥teus 
sreCst tite LB whee ous a8 ,eNViteout MTS 4SSeSEOUS 
oeeCit om Cs .elotreds 8B ,d1VIE CUS avots ,deeerees 
SMeL tt me Ca .@esSteu! aS 2 seh o0u MUU ©, ¢2S%E 005 
peru ge fo CB NEGF Ong BS ywiedtous wWte eo Peepp woh 
Serine ote LB ,paveteuy ae gfivatecu VB peenwreus 
yweoule ¢ CB guerre bare Lc ey O02 Oe erude ty 

shoulder oe Le ebbUNe ery ae yfulereys eels e,nanegeue 
8 ee (8 ~eelateny ae ett steus oulls =, fesiteua 
veel lt oa CS obtedpens as ,vehatetios avete Oy 


man 
vis 


arrcit 
Soecit 
SHeCit 
sPecit 
wecit 
S¥EC It 
Sle 
vet it 
serie 
srecit 


as 
vis 


Svecit 
SveCle 
acc be 
svtche 
weet 
shecit 
SPECIE 
SPEGIE 
seecit 
Srecie 


mas 
Pta 


Nhec it 
SPoule 
SPEC It 
SPECIE 
SPEC IS 
oven le 
Seecde 
specie 
srecae 
SPECIE 


“oe 
vis 


sPecie 
See te 
SY te 
orecit 
See be 
SPECIE 
Svecie 
grecte 
sveecie 
aeeule 


“as 
els 


weas Pa ee 
visre elevbre ad 
ries ebbsare f 


ev a8 bt .O2o0tees 
o2030C003 
O51 16008 


rT) XG 2,400 St 064 
ob Ptef 004 
oF eleeul 


to 


wy 


ue 


ed AB teelteeout 
olviuteul 
oPS20E C08 


tw 
tue 
n 
ne 
tee 
neu 
'¢ 
uv 
ue 
a 


és at Py navar out 
ehbircedl 
PLDs) 


te 
te 
e 
ne 
Un 
neu 
re 
rd 

a 


es AS Pebusveeud 
4as ge uveeny 


VIS ,Viasteony 
SPEC KE ues 
peeCie le 
setcih or 
an ard 
ld oS 
SPREE omer 
SPrulh one 
dseeCit i 
swoelh ove 
Sercit ® 


> aS te/Suleevt 
ohVnNe ent 
et¥odeeny 


'tctals 
ri teags 
Arye 


ms LSet Te ees 
oFOditeoe we 0.22068 003 


dEL te 


wnve 


ts 
te 
cs 
te 
ce 
te 
cs 
ts 
ts 
ts 


20108003 


oPOSOE 200 
eb Ht oul 
ob test eos 
ofet2¢°9! 
2200€ 063 
Levet coe 
oS 200F 400 
otagseoos 
el UHRE ONG 
0 1800t 965 


We 8.582% ¢0t 


OeLTE 


ams 


ts 
is 
Cs 
us 
Cs 
Cs 
te 
Ce 
Cs 
ts 


Onis 0,253 001 


Sue ,06seh es 


1s = .2te3ECes 
XG, S0K0E CES 
8B datogoee 
Rs .240bkeee 
As ,d00eGoes 
Me .2toigeee 
uBR .29°SF 000 
RB LL Ieote0S 
as ,vleetove 
se Daeetete 


Puts 3,9396E°0! 


ehOOSE 000 Ha @, PASSECLS 


odtesteus 


o2¥ett oud 
ohSUTF CU) 
oteotteus 
024596 201 
efOS0t ens 
01909St 000 
242008000 
016556005 
eb 1SS¢0uo 
el 00F ONS 


ke 1.0858t0} 
eOV0AECOU HS ©, 1000te0S 


beLvs 
anus 


Cz 
ts 
Cs 
ts 
ts 
Cz 
ts 
Ce 
Cs 
Cts 


2p see0s 


ee¥o se eta 
eOS7et oO} 
01503003 
2€OSGE oul 
2 355Gb 003 
eS ¥uer oto 
2006 Oty 
e2OSaR00§ 
rei riss ty 
ef tOvteus 


We 1, FUNGE oes 
eo UUSEOUU ME OL 1OeEous 


URLYS 
anus 


Cs 
Ce 
us 
Us 
ts 
Cs 
is 
Cs 
Ls 
ts 


edi sereus 


Py ah eee 
ed IV2F oul 
el errgeus 
oeaTurent 
o3teSt 05 
oS 8O3t 000 
eM PONE oUt 
2 COWES CUS 
ehh Leh e's 
ef beckons 


me LPO 
eOPUTE ONO MB ©, PREM ENS 


OPLVS 
ahs 


\s 
12 
Cs 
(s 
le 
La 
ts 
(3 
Cz 
(3 


ee Seteous 


of SSet 06 
ebednredt 
eV STSF end 
eC URS ont 
ah sabes 
Peay 
oHeG(e oud 
od brtteny 
eb i¥erens 
aliutrerdcs 


es t.osusterd 
ehhehterd me e23iseeus 


ELS 
megs 


eeP$sureous 


B28 ,VSViL CC? 


aa .ctloteor 
BS ,2vo2tens 
BB geusvEeve 
BS 4 4458E 000 
ae SUSIE oeS 
RB dtd3Eoee 


is ,2veeteeo 
Be ,200eE obs 
1B 4 78OSte0e 


ae ,5aN0toue 
Pale «,5765teu) 


Sxe , 95500001 


aa ,2tigeeue 
a8 = ,24983€°0) 
Re ,209@oo2 
RB 4 28S4E 400 
Rs .G@171E°0$ 
Be .ehustevu 
as ,28N5te00 
Re ,Se8ateos 
2B 2423805 
xe ,S888E oe 


Pris a. SA2eten} 


San .VSieboul 


Bs .2102E 000 
ae, SuMteony 
ae gendleoue 
as ,duedtoud 
Ne, $000G 003 
ke ,20N0E 000 
Awe 4 29@RE 600 
ae Sb Feeus 
1s gLeSME CNS 
as ,SanBeeve 


Pris e,eseoteou! 


Sus ,VerTeedy 


eR ,2107TE CO” 
as, S033E ou) 
18 ,2laateve 
QS pdunut oon 
ta a Sohueonds 
ae ,2luakeuy 
AB d9RGE C00 
ts ,2aVeto0S 
1S LLIESF ous 
as DAVE oUt 


rHLB 4, eddeteus 


Sua ,@25¢teul 


te 2.200theo3 «= a 2. Bnven ces 
ats _3580teo3 tance 6, 
Sumputs .e3s2ECes 


OPATs ©. 21638 oo 
edits § 28eseete 
avuts ,S008foee 
MIs ©, 460%ohs 
ROOTS <,S000Eos 
s00ts 
avOTs 0, 
abuts 


Ta 2.29006 003 Pe 2, 9415eeu0 
wis ,$023F003 8 thas 8. 
Sushuis 47996405 


eurta ,27sitees 
21d Os 
o. BH 2tooe 

-*5eeteoe 


or “cotous 


eiite 
ROUTES ©, 8595E0R8 
abuts @, 


HOOTS 1 7ECECOT 
ALTE ©, 7000009 
AvUTE Uv, 


ts deadditees Pa 2, /400t eee 
asa .870etees)§«=6 Tavs 0, 
SumvUts § .S147ES0S 


HLUTS ©0608-2905 
BLUIS .0808e0S 
BOUTS 20038903 
BOUTS ©) 9888003 
AOUTS ©, 2020Fe0§ 
abuts .19Seb~oo3 
aunts 0, 

AUUTS ©,9932E00? 
POUTS .Se0NEes 
wO0ts 0, 


tn ¢,25128¢03 Pe 2.5080 +00 
wis .3seeteus§ «= Tar: 0, 
Sumputs ,5905t005 


AUNTS ©,23016°05 
AUUEB oh SO1E COS 
nutes ,e0eokens 
Autts ©,2053003 
ADUTS ©.26 7003 
avvTs .2098be03 
avuTs oO, 

nVUTS ©,ub)1£ed® 
SVUTS ©. 310M COR 


aluts », 
Ts 2,20616005 Va 2,5833e000 
wtz =, 3AsAE CO} Taws 0. 
Suvcats ,abvoteds 


HUUTB ©,4 490005 
BUTS ,ELOOEogS 
aunts ,820sbed§ 
AVUER eb hVGbous 
suerte w,heiecods 
AUIS 612 SoOS 
euvuls 0, 

aula e275 7Es0e 
SUnTS ©, LQ0TESOR 
ruts vs 


te 2.1Sedsteous Pa e,ovlateve 
wta ,295hFo0$ Tanne ©, 20772005 
SumpuTs ,eSeséeus 


AALL TEMPEMATURES 1000, KELVIN 


AUP aMT Cath PARAMETERS 


Us 2.1nettees 
eee, 


UO 2. 1bbek eos 
ve @, 


va 2.189% 00% 
Gee, 


UB 2. 1F2ebvE4 
us 0, 


be 2. 1000€¢6s 
Qs 0, 


Va 2o09S7E 008 
wa o,10s0t ees 


C0600 000080 a BOOS O Ee OAC R ATE EES SOO OE VOTE OEEOSTOVESESS FOOT TOS TY SET OE ESOS OSOUOSES SESE SSOEDESSSE EE SESHOSHESSHSSHSH SO OUREODRSOERO 


adi hat 640 im€lAte = .elesb ood Gris = .3Fust-0e «Ube | JL Sae Oe 
ee ORE ate inf (ACB Pre Ld re | Chaos o4Oeebed? CC sCK IS oH AORE OMY 
este papers fhuveBeus s.ee7 Fees 


Iss] 


Se OD RE PEACE IEEE UIE ISTRY ARINY Si 


on ieee 


q 


sted kite Sandie tr malate hal oe RAE: 


. ~ ame te em te et 
ae 2 cca = qe eeinensnanemane emanate merase irae uantenr ee — ane tenn nn meena Ne cite oe 


BUUNDAKY LAVEW PRUutILEs 


& 4 ust wit usub 1a vee ' 
eS ofeverees Geter sun 50aSECe0 e2ivetoes 
va i o2#00toos 087736900 o@0oTteey paleteess 
a i 0 30008908 0 €855¢ 400 oVeeatete eO276to85' 
eS i oS0U0E03 o89a3ECus oleosteos ob S24Eo0$ i 
RS eS00vt~ous oSusut egy eluOng ous el @Saboud ? 
f eoueroeed) oSt276 000 ob /AVtCNl ot2snten? * 
j ! o?000b°u) oS]22steuy o2070b 003 01430tou2 
ie i o800UEC0) oSS17E000 odSeSieds oloseteos2 2 
5 ! o¥OU0Eeo) 0541 3E +00 22009 008 ol8Sdtou? 4 
g | S huove eno 35$10F 000 s20S0E ous s20a7t a0? : 
FS ol Suut oun 259%et O00 oevetoul o@Vletoud - 
‘ ofVOhESON obe7at oud oda 3M 003 0 Sant ove : 
ar o€S0¥E 000 o0%Set ou 200150008 090 35t002 ‘ 
% o3P00ECU0 e7S75t000 ot 00Eoe! edt tetoe? 
oe o3500R 000 7 788€000 o202E 008 eBiedtood 
y 2@000E 000 81726000 of0eSE 002 ot0S0Eo0s 
BA eS$0ut euU oBS7SEe00 oS 15996 oud otdeStous 
4 edvovecvy o8RaSt eon alélveu2 eles stot 
q eds0UbeuY 293 30Ee09 ol Savten?d oft iteovi 
* 060008400 oO38uteu0 olapetee2 0 S527F Ot 
t @S00E 000 eVS9GE CUO oIS@rEved 094050018 
4 ; 27000E 000 29769 900 ol2 008008 o@2aat oot 
ae : oe TSuys ove o¥Oule eon oheudst ode etesatoue 4 
@ oMQ0NF Suu etOoteeul ol GuSeeue oS 9u2E e008 
; oF SUOE COU o1000t +08 ofOt Seoul odhgigoge ; 
ie 29000E +00 o1008E 00! oh i TEsoe 4. eee ! 
x 299006 400 of O00E 003 o&220to02 esototees a 
x of Q00RCUS of 000E C0! o2S24E 002 ebtuetent i 
3 ¥ 
E. i Seecae cu CE 28antoood «= FR CgeN ORE OUND «= ALU EB ©, 20078 60% : 
B. Breuer ute CB 4ORGUE CLS ae y MUogeul avuls ,1007L00$ 
BS i ovtcle or Ce .WeSer ous a8 .leaagoue evuls .7eeTeoee 
ee) SPECIe m2 Cs ,aeloteul BS 42ul2ee00 edut 
e z arecit um be 917418903 Be .2306E00S) «= UTS, t 
A SPECIE met: CB g1BOSt COG We .2300e¢0u AvUTS ,POTOE Oe { 
c svecie fe Ce .a206 Oey Is ,2vOSeeu0 suule uo, i 
" Swecde ou Cs a weuureous RBC ENSut ous wbuls eh173tove « 
: sveutr ove Ce le btene as ,a512te0o BUTE ©,43541t 008 4 
‘ = srecle oa CB ghidurets Bs 4 244eb oun eyure u, 4 
TeooS PME APL Its t ‘8 Be grwns 


oe LUCML wad HmiiPrmt Les 


UG SeUIUSeVS VS rewhnbens FEKRSTIEVE FUP ena WETMILYS Wee MeVEEEOAD 128 BL AP24CCOl Nee BH 2Peteo) 
wows MARTILLE PHASE HaHERTICY coves 


POWFTELLE weuur ’ e 5 « BH 

VUemSimbae veLueir LeoVukeudS 1eebirerS 1,930600S teazieeuds i 

4 CHYSVOSTMEAM VELUICHTY oS,Po0keOt WP 9S9R 402 04,985 e02 o¥,alSEoud a 

4 PORT TLLE ARYL de TH, 1.e2nt eon hentbeeut Te453t 00) = SPSL CUT q 

: Pahl CLE Heewe mating Ser2MOOS «Porshe ee TISF ONS «=e PNR ENS : 

’ MAb TEele veudtie 1,2uat ous TevISFolS GL hIIb ONS S.2esteus i 
VAWTILLE ome cium btue €.eSvt eve ‘ 

i woe THEAMLINGS é@ $s 10,2090 4 


we CUCAL GAS PHumewdIESs 

UB AMUPEGUS Te PUSH OES CF oOLVE S,CASE UR WEYETLOS wU,8 ©, 90PF OP] Ves toourFedu w2S 1nazéend 
weowe FANT ILLE Hoadh bauPERITEN coece 

PARTILLE Gree ' e s * 

Vuanalegam whirnetty 1,092eeud Ieohieeud F.598t40S F,02sb oud 

CHUB—eST HR aH VELILETY es ,bUSECUE OS.SUE OU? eA,wttFeued 91,090b 008 

MORTECLE WH rtemery atte L.2Gueus o,FurFeuvd 4,44 401 «2, 3neteut 


fhe 


eee 


4 PARTACLE Te rtmatuwh — abSestous  eeerrteod  arhubens ey TAME Oud 
F Mabie de Ferearts bedulheous be lvebeuS Veev0Eeuy eetidtots 
POMUILER eeeee Tee bead eodeut ony 
soostebLaal burs 4 ys Wess 


eo LUCAL GaS Karivg elles 

UB CeULIteUs 18 2,u5UPOGS VENSIINS J,0AZE COG WETAOLUS W.8 O.d7?to21 V2e 2,a61FeuN N28 2,050b 000 
@ooce PARTICLE HONE PRILMERTEES cocwe 

HATE le Ueub ! é 4 « 

MulanStmtay vecuctte Te@¥EE CUS Lebmlt ous 1.S54t on§ tedeot ens 

GuddedINE An VELUCTIY ee SYNE HUA OBOINEOUS S8.eddb 003 oh SedteUS 

PaMVLELe Me TAML OS aU, Te@33E000 Ow, fiPEeud 1,487800) =o JesE eo) 


q Pawlicre Meawewatung eSPSECUS 2.egbEsus 2,/108003 2, 7066 603 
b Paetdcde Ubwsite 1.P2ebeous OoMISEOLS «9, 10SEMUS =e FURFeUS 
PAT ELE made ceees bbe Pedant eu 

wooS ink Arn hues v 98 lusdvtus 


om LIC Oe GAD Hmibentgto 


US PeMtUEOUyS Is 24540003 CENSET TS S.OPSEeUN REYNOLDS NU,8 eBSIECCL Yes 3,280Ee0U NPE $2766 04 
woowe HAWIICls wrice whiWhnT IES eavce 


Parrderl Gand 1 ¢ 5 « 

DuerStwe sr wth cite PeMwlb ees Pannibenys Ledsebeus Tesenteus 
Cor sdeniant ® wthieddlt sazmagt ers etycurber eof feuteus ee ntaters 
PaoVLLe ME eet > te VadSut eur tone eny beGudt oud ee 55h Revi 
PAWEECLE Ieere nadine Poet org 2.02903 «SPITE OHS «=k PBHE CUS 
SawTILLe ve eotis L.2eutous = oebevewuS 9, 1MuRO0S =e MBBE COD 


POWEICLE Meme ise bhi Pele oun 


thems were dail 


a er eo ere 


“owoS int art fut 8 Ss ss tu, eevos 


wo LUCOL GAS enuventies 


WE S.9LILOVS 16 Z.aDeR OES HENSLIVE S,VOdECHR METMRUDS MUce SeHSEZE TZE C,CUSECED Ze B.000t Ud 
eooee PARTICLE PHASE PHUPEWTIEY seoce 


DANTICLE wtuur ' 

Ovahstaram weuocsis $ MOTE O05 
EMNSRedietame vere her oe, (NS eu? 
PawbiCre meres ai, fedsdtevy 
PanTICLe ie-Penateat eod2St ous 


é 3 
beemarery §=tS34EouS 
*eolleeus ef, 8aSb 005 
MH meat ent te sutted) 
Cooder ons 2eTidteus 


Pamticce venStey UeeeIH Oy = AINE PES «=—«- BW BADE COS 
POUTICRE Heme EUS FLOR 2.939t 006 
ooo tie AML (nts o 38 10, 20004 


ee LICOL Gad Peurent les 


“Ul PeVUME EUS FS Poa MR OHS Cb eSTENE QISTReOe HETHLUS 
owwe HAMIILER PrHAne Farben Its evoee 


PANTILLE wavur ' § 
OASIHEAM weLIMdtY beret ous ,ceneteuvS t.599005 
CALSSOSINE AM vere’. eA,gidtenr?  -P,1abovds oP AuSE Ong 
PARTBLLe HEV DD VetPat ers weOUvAF enn toSsnr out 
Pawttcat Ve rematinn Perdtbor, erdotevs eltarevs 
Mantle cevodtr terete th HeMHareny VeAuiteuy 
PAmIeLe caret tus uve e Shee eee 

woostnrartines a 98 1e, 405 


Om LIL AL aS Pormemdles 


Us 2.UNMeeS F8 P,edseters It vn, 
enwoe ranlitit rhadt erurtHl les eevee 


PAMTICLE beiiuy r) 

OVerasingam vearCity V,H9CR ONS 
CHUSSOSPrEAH VELILEIY ef oeteens 
PAMTICLE HEV eto tr, Pecdht eae 


Soladteons we reins 


é 4 
LeehrFers 13 CG 
2. UOVECUS OD,07MR CES 
pew ery Fovstbevt 


PawhELLe Te -Heneieme SeDAVe Oud Comerheudt Pelultous 
MONTILLE Urnolie VekSHbons — ems eS WP tHLons 
PARTICLE mumtiditm been @edeit ory 

ooking ary ties 4 98 10,2040 


@e LUCAL &43 waurentses 


US P,VITE Ors Fs Poeutt Ord UbndtiVR VeCPsreus wEVPULUn 
ewone PSNTILLE Hreadt Pel HtelItS esses 


PATEL Ee Unter ' 
Quart stekaw vetrt tis VeOYit ons 
CHUSSeSimgae wert bit eF,eest ons 


e s 
Pemereoud | 1,964,005 
me HOSE 008 


MONTILLE Bt Heros has, ys Ce Ph S| 
PORTELER FE rer tant Pedtht eos « eemcibevd  seFvoteus 
Pawlarce wae ager VeAotbenys eF.95rterS Beuprrenys 
wo Pantene a ee rhe cat Tb ony 
amos tie AUR ENE Y 98 1G,ev90% 


@o Lib OD wmihewlics 


UB 2.UESEOUS FB A ,wadt ors LENSERVS ¢,00mteus RETAILS 
wooon VARTICLE btAXL Ome T EES conce 


PArTILLe Umenim i) 

War >INb As velLucats 1omdee ony 
CMUSdeSimeae veLiCtit et oleceus 
PAWTILLE METAS wh, Soeeit ean 
PARTILLE Te mem atime Sobt Steud 


¢ $ 
dehrnrevs = PedSat oud 
HH,8Coe ous eo. Seseens 
@oSISEsLU FeUIAt oul 
P.Piorers ee Puakons 


PARTICLE vet vide Peddetous heWhavors 9, carers 
MamlbLee ter bee ture Per Te ony 
@enoystnranry | ot 8 Tiny Torr hh 


Oo LULAL FAD Payrtadicds 


UB CoVdTeeUds Us eoehHeers FRISAIVE 2.HeTbeva Me TALUS 
coeew VONNILLE raat MaibertIES coves 


ParVELad teary ' 

VUArotar ar eburc ate Lovers 
Ca Sregtnrd eee hss ee tuteons 
POOTILLE we eeLUD we Veeder eun 


? 4 
bemvinerS Tereor ous 
ed PUA Oe ON frre ens 
wo dadse enn ba Ste eut 


PARTILEE Ve hema Tiint PeSUG CUS «=P owEEE CHD «= Panvnt eas 
POmTTELe Uewdptr PoPStie ous PULP ooh) @.$eet ony 
MPawt year sete few Flee eaPtSt one 

wee Name Ars gies Yd of te ayusnd 


eo Linde Wad om rcmdiad 


UR Pavel euS BB eWsmrons MenSel es PC avPleoue REYAHLOD 
weemws FAUTILIR PRASE Priurer tli S escoee 


Pam TLC Ld tomer 4 

MredvIngas vty Gily Pem@re e005 
CWrnventwe ar vee Gaels ee lebarons 
Cy ee © 2 ok VA hhh enn 


é 4 

1 hVnt Oud Vedeut ug 
oN. dobet§ om ete oul 
me imi bev t, WMat out 


PAmE ieee Les mea tung Perdleeus eemiah ous = A,puveeus 
MORN GCLE Ce oSEte bePShe ers TeYosredy 9, 3NStouS 
Parl iL he sehetedsee Ohad eH SSE oun 

eoostapa gl ve Ve Sa Legeraegs 


we yp tak ho Mortem ites 


HR AecANE OVD Fd eowbei ands ct SENN Germ ttona Merrit oD 
wowoe vamiiLit Yrest erubent IES cacce 


POETIC LE Umiive 1’ 

Vieroiagar weLoudly Lebar eus 
Carssestega ves bate eb OSpbeus 
Panter abd rq ON ty PP hud one 


é 5 

Pemvot ony te SSet ous 
PheIWNE OS OP ULE Ons 
eo SVPe one 4, done oul 


7 Pee ee Oe ee ee eT Cedtst oud vemtiters Pemvetend 
weelicae ceoaale Leematens Febetbeous 4,519 ens 
vemlitig oo te eh te P_AMeh sen 


a 
1.edot 609 
oF STE eos 
ea Sout out 
ep tastees 
2,88 205 


NteS B,THSH OP) VAS G.VELCOEH HPS &,POnkCEy 


a 
1 S2SE C¥y 
oe, MIE ous 
ea Sant ont 
celhigons 
So hFEE HOD 


NueS Gefdeteet VRE SPSS CUe Hes 4, Tout coy 


a 
3 .02Stous 
oS, en0teed 
deSSsTeeut 
eoTTeeeud 
eas ous 


WieS MeT20R021 128 HSSIL OOH WPS OSIKE ONO 


a 
1,453 009 
oR 1OIEons 
2, S03 oot 
eo Teitous 
eV Ar ony 


MUS Bcb70F C71 VSS F,SVOECOR mee 2.9198 000 


‘ 
Veasut ous 
oa, beltous 
2. sMokeur 
eo FPMevs 
F4ASE OHS 


WtreS Sebsotedd VER VesVpbedy Mds O,1atk eve 


a 
Peasitens 
PO mieceny 
ea Sudeeud 
? Tooteos 
eo ¥Sereds 


M.S Sod¥redt Wes OOtatevse KPa B.e71to1n 


s 
Peed ons 
ee,osefeut 
PeANht oul 
ea Taoteus 
2 Vengeas 


NHS HNIID OPE VES Vmaukehy wes Veet enn 


“ 
bee hae ous 
es SVat cus 
Peedi oud 
Polite eud 
ee VINE ny 


B-13 


oa Le tae ae a ll 


- ateu cp Laney Ser ET, | St tn ane nav nae era mms Te a uo by RE 
poemewr AS IGS EMTEZIORE SSRN STATS TS, ok eee a ae - 


\ 


wood tht ant int 8 1s «$s tvedb¥09 
Oe LULAL BAD waieHetl its 


US SUTEEOUS 18 SeaAVELS LEUSLINE 2.H235EoUs HEYRULUS RUS O.eSeteds Tee 1.C00F OOS Nes 1.ONtF O08 
evcce PARTICLE URAdt FHUPERTIES cocce 


ehh 


% 

cal 

PARTICLE GmuuM a e 5 a S 

OGenstteat veLrcity $.900bees 1, 70SECOS «= SSCECOS «=. 4 ER ORS, ;] 

CRUSSRIAL AH VELNCITY oS.SNUECUS 0.15008 98,8428 083 98,0008 405 * 

PARTICLE HE VAIS Mts, VePSTECCO = = AOE COO PeSO7TECO) §2.788K 008 4 
PANTICLE Tht watt eedO2L0US =e eICESUS «2,09RE CES «=e FOVE CUS 

PantiCee versity VePodtoOs 7,163t005 VY e2atouds §,008t~es 4 

PAMEICLE mumentUM FLUuA 227186406 3 

ewes tae Am ines le as 10,20005 Hl 

oo LUCAL GAS PourtHtits y 

Us 2.000R US Ie Poueseeus DENSLTYS 2,.75ifeva REYNTLDS NUS O.Selfedl ves t.8S2Foes Nee t.1dth CO) 4 

wowce PARTICLE PHdyt PAUTENTTES coece 4 

F PantiCae exuur t é $ 4 4 
3 OuenStwEar verulaty VePAIECOS s.JomEeuS 8.599E00S 1,081teNd 


CHUSSOSINEAN VELOCIIV o4,S378E005 98,0SAF OCS oO, 4OSECUS 91,0008 008 


8 
F POWTICLE WEVNULUS WU. e2StECOU b.2uSECOO §«=tSSHECU «= 22S AE Cot 5 
is PANTICLE ToMPENAI UNE 2494008 =e HUGE ONS §=—2,wOIECUS «=. 2. 788 003 a 
a MAWTICLE VewSity L,Plateus T.etetevs @feetoud §$,8048005 a 
a PARTILLE stimtvius PLUS Felust oo f 
7 oooS int AM ines 1s os 1o.209us 4 
. we LOCAL LAs PrUPeHTIES 5 
2 Ue 2.O5eteuys fe g.888EC03 CLENSITIYS 2,000teds RETAULOS 40,8 O.270te2s V20 1.237600) wee 1.291860) 
be eooce PANTICLE PHASE PHUPENTIES coece 4 
Ss PARTICLE bMuUY 1 e 3 r) j 
% OUMnstWEAr FELICITY te Vare ons e716 ONS 1.507teUS =. 1 48M CUS : 
aS CHUSSeSIHeAm VELUCIIY @$,.5006003 O,V0FE COS o1,80SE CUS oh SO7E CUS 4 
a PARTICLE MEYNILUS mu, 1e@9SEC00 0,2826000 f,32eboul 2,151tee! 4 
: PanTICLeE HeeveRatunt PoAV2ZEOUS 2.0U3t903 2, 090F 003 = 2, POOE OS 4 
‘3 PARTICLE venstiY L.PTTECVS =T.S3OFSOS = 8 00SEWOH = 3, PASE COS 2 
a PARTICLE Plat sium reUK eoMSSE OCU | 
a woeStutact lies te sa 10,2094 ‘i 
By wo LULAL Gad PaurPenties 
ae Feeuy fe eg SVGEOUS LENDEIVS 2,D0BEC0G WEYRULDS NOS O.8O2EO21 V20 L.32UECOs M2B 1.201800) 4 
eooee PANTILLE PHASL PHUPENTIES cocee 7 
j © PaWtICLe Gaui t e 4 rt) i 
VURNSTabde vELUCite 1. Yemt ony deFete ous 1.572€ ous 1,ente eos 7 
rF CoysSeygtecam vecilstr oF uasteus of ce( Foca of, e9dt ous ef endcous \ 
q PARTICLE MEVAHUd 1, FeSSTEOUY OeeHsEOUL 1.523800) dd LOSECUL 4 
“ OAnTICLe Temwewatune 2.47% VS 2.582003 2.072EC0S 2,7426 003 ' 
PanlICce Obwssiv 1,274¢ evs 7.d358E 05 t.09%Eoue 3, 5038905 * 
POMTICLE MeeeNtum FL Seong out 4 
c. SIME AOL ines wv Ss 10.20905 
; »7 LUCAL UAd PmMEW Lita 4 
x US 2.07 eUS 18 SaSeee HVS LENTEIVSE C,ULUE OH KEVAULUS NO,8 OoddSted? Ves 1,050E00) Wee 1.373600) ¥ 
wooo PANTICLE OMASE PHUPENTIES cove i 
; PARTICLE Gruuy ' e 5 « f 
4 OUmasiee ar wersibads 1.9$me 905 SeT3eFeuS = 169808005 = 1, WP SE OUS ° 
4 COUSHeSImEOs VELOLETY em nSures of, SOsbeia e2,0UTEOUS og, andeoiu 4 
PARTILLE WR dweES we  Tedfoteun oodsikone 1,5ESROUs e,deMb oul i 
q PAWTICLE lemre watune Coto eos adhuteuds eeolereds ep festous { 
Panttrrat vevwstiy $eeVeto0$ = B Jantebs We2ssEous a, 7eAcous 
q PARTICLE MAF ele Boe $e5U Se 900 4 
¥ eoosihhanints 18 98 10520905 
ae we iCal bas Pomentien d 
' : AoASHovS FS eSdieers LeUdPes VS, 52it ous KETAKULUS NU,8 @,uUSSEOeL F2E 1,50Gk001 hes |, 0S3F Ov) 
wocee HOWTICLE wtaSe HAUMLRIIES evoce H 
PanticLe wave 1 2 3 e 4 
DUAnASTREAM VELULEEY 1,99%eeu5 0, 0. o. 
CMLIRTEOT WE Am VELILETY om igubens 0, o. o, 4 
PARTICEE at vratns hus | LeOTSE OOD ou, v v i | 
‘ PAMTICLe He beaatout eousseeus 0 ve 2 
PARTICLE Ubuolte 1, $aapeous Ve Me Ve 
“4 PAWEL ER ctrme bee ete Pedant oul 
ry eoodiht at hues iy ¥e 90 e2u9ys 
i we LAL ad Pw thew dats 4 
Pelareens TS eoLVIL ES PRNSPUTE Ce 1MbP es MEMOS M1.8 S.HHIECP] Ves 1,580 00) wes DSPat oul . 
ewowe HAMTELL? Pract MRUPCRITES exces i 
.. VaMVICKE wMnr 1 e $ e j 
Ouamdslnt ae veuoclir Ue 0, 0. o% ? 
‘ CWNQdeSTabam VeQucIty 0, 0, 0 o. , 
b Pawhtere Meee tE eS wer, he Ge Ue Ve 
. WAWTILER FR ob wating o, We Ve Ue 
i Maelicee Gtr gate u v, Ve a, 
A Paniiece sume stue bLua Ue 
eooSing AMLINES ev $s 10,20405 
we LICAL Lad PmuPeHTIES 
U8 Pehulhous Fa P,Pvibevds CE«eSLEYS S.1eStend REYNILUS M,8 Smite eed Ves Vemerl Ab] Mee 1,SASt oO) 
eowoe PAnIILig CHOSE Hmiimentits eoove 
4 hantlice UN ' e $ “ 
4 Dummatweae viculite Ue O. a. ve 
Ceusdsestwtat veructit vy v. 0. Ve 
PART htLh ME rks HU, Ue 0. Ve oO. 
Paetidee Meee wating ve 6 v. te 
Bar tICee Uevotir on We *, Ve 
Cae ee eo oY fe 


BULNIART Ob wari leit wedse al Ss PVOECUG P.H7SeOOL DesbIerud I Sebeus 1, feet eos 
BUUALAWT UF bom TiCLe wmacst al Se $1,0109O f,uGobedd 2,GSdE oUt 1, SeGeood 1.553 008 
BUUNLANT Ub PaNTICLE HmASS 41 SB 14,0096 ESTSECOL §,0G7EeU §,.4088001 Jb, Sent oo! 
BitnePARY ob wanTiCes eadst af S8 1,67GHH T.Shneout Leemetout towbPeeovs 1, 3Peb on) 
BUUAYARE OF ran Tit Ee vaSh Od oe ie .@Sid0 PSateut bedPoreud Teehstenr te sonteous 
Poteet ue hawt ke © doe al Ss le miadm F,SHur ory Feaeseedy b.4eSreul 1, eust 608 


APPENDIX C 


Z 
: 
Z 
a 
By 
> 
& 
¢ 
2 
a, 
& 
=) 
5 
: 


(is! 


Te ee ae < ae ee 


APPENDIX C 


SAMPLE OUTPUT - EXPLANATIONS 


The first page of output contains all the input data on Cards 1-5 and 7, 
followed by the species identification data (Card 9) and the distributions of 
properties along the initial orthogonal surface. The units of each variable 
are identical with those used for the input data. These are followed by the 
axis location and the nozzle wall contour. 


The standard printed output at each value of KP (Card 5, Cols. 11-15) 
or DXLSS (Card 3, Cols. 11-22) gives properties for each streamtube, from 
K (streamtube index) = 2 to K = KMAX. In addition to the usual output* i.e. 
X, R, PHI, T, P, U, the following properties are printed: 


MA Mach number 
DELY streamtube width (cr) 
H enthalpy (cal/g) 
HT stagnation enthalpy (cal/g) 
TAW shear stress at tube interface (dynes/cm*) 
Q heat flux at tube interface (cal/cm*-sec) 
PT dynamic pressure (atm) 
RHO density (g/cm!) 
SX distance along streamtube (cm) 
SUMDOT total mass flow bounded by streamtube (g/sec) 
Cc species mass traction 
Xx species mole fraction 
WwDOT 2#pecies production rate (g/cm? sec) 
only for} Fl total mass flow of species up to present streamtube 
viscous (p/sec) 
flows 2J species flux at tube interface 


* Xand R refer to the coordinates of the outer boundary of streamtube 
while the flow properties ire average values across the streamtube. 


q 
4 


Boundary layer parameters, if computed, are printed between the tube 
properties and the composition data in the last streamtube (at downstream print 
stations). The following parameters are printed: 


RES Reyncids number based on distance along wall 

DELTAI boundary layer thickness (incompressible) 

THETAI momentum thickness (incompressible) ‘ 

CFI skin friction coefficient (incompressible) 

UFI friction velocity (incompressible) 

DISP displacement thickness 

DELTAC boundary layer thickness 

THETAC momentum thickness 

CFC skin friction coefficient 

Hi2 shape factor (= DISP/THETAC) 

QWALL heat flux at wall (cal/cm*-sec) 

TAUWALL shear stress at wall (dynes/cm*) 

U/UE B.L. velocity profile (UE © velocity in last tube) 

T/TE B.L. temperature profile (TE = temperature in last 
tube) 

Y/DEL location in B.L. (DEL = DELTAC) 


If IFART = 1, particle properties are printed; first, a NAMELIST of 
the input data on Cards 15 through 19 (and 20 through 23 for NC = 0), followed 
by the initial particle properties in each streamtube, The standard particle 
print gives the downstream and cross-stream velocity, temperature, Reynolds 
number and particle cloud density for each particle group and total particle 
raomentum flux, for each streamtube. Since, tyvically, the limiting particle 
streamlines are within streamtube KMAX the outer tubes will contain no parti- 
cles lLimliag particle streamlines are noted as, ‘BOUNDARY OF PARTICLE 
PHASE AT S$ =X. XXXXX', followed by the radial position of the limiting parti- 
cle streamline for ¢ach particle group. If the particles are imtially in the 
liquid whace the program will print where they stezt to volidity. 
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